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Abstract. The theoretical predictions for the subtraction constants lead to a very accurate dispersive 
representation of the tttt scattering amplitude below 0.8 GeV. The extension of this representation up to 
the maximum energy of validity of the Roy equations (1.15 GeV) requires a more precise input at high 
energies. In this paper we determine the trajectories and residues of the leading Regge contributions to 
the TTTT amplitude (Pomeron, / and p), using factorization, phenomenological parametrizations of the nN 
and A'^A'^ total cross sections at high energy and a set of sum rules which connect the high and low energy 
properties of tttt scattering. We find that nonleading Regge terms are necessary in order to achieve a smooth 
transition from the partial waves to the Regge representation at or below 2 GeV. We obtain thus a Regge 
representation consistent both with the experimental information at high energies and the Roy equations 
for the partial waves with I < A. The uncertainties in our result for the Regge parameters are sizable but 
in the solutions of the Roy equations, these only manifest themselves above KK threshold. 
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Low energy pion physics has become a precision labo- 
ratory: the chiral symmetry properties of the Standard 
Model can now be compared with low energy precision ex- 
periments, using Chiral Perturbation Theory (xPT), dis- 
persion theory and numerical simulations of QCD on a 
lattice. In principle, as illustrated by the prediction for 
the magnetic moment of the muon, physics beyond the 
Standard Model can show up at low energies, provided 
the quantity of interest can not only be measured accu- 
rately, but can also be calculated to sufficient precision. In 
this context, the interaction among the pions and in par- 
ticular, the elastic tttt scattering amplitude play a crucial 
role. 

The Roy equations [Ij provide a suitable framework 
for the low energy analysis of the tttt scattering ampli- 
tude, as they fully incorporate the basic properties that 
follow from analyticity, unitarity and crossing symmetry. 
They express the real parts of the partial waves as inte- 
grals over the imaginary parts that extend over all ener- 
gies. The high energy behaviour of the imaginary parts 
thus enters the analysis, even if the Roy equations are 
evaluated only at low energy. More specifically, the con- 
tributions from the high energy region are contained in the 
so-called driving terms. Near the tttt threshold, these terms 
are very small. Accordingly, the results for the threshold 
parameters are not sensitive to the high energy contribu- 
tions. The driving terms, however, grow with the energy. 
In [2111] , where the threshold parameters and the coupling 
constants of the effective chiral SU(2)xSU(2) Lagrangian 



relevant for tttt scattering are determined to high accu- 
racy, the Roy equations are solved only below 0.8 GeV. 
On general grounds [1] these equations are valid in the 
regiorQ s < 68 M^, i.e. for y/s < 1.15 GeV, but in order 
to use them we need an accurate evaluation of the driving 
terms. This calls for a better understanding of the imag- 
inary parts at high energies and motivates the present 
work. 

For TTTT scattering, several sets of data at high ener- 
gies are available from indirect experiments [71113). These 
were analyzed in the late seventies, both in the frame 
of the Lovelace-Shapiro- Veneziano (LSV) model and 
by means of various dispersive sum rules ^15,-22j . While 
the determination of the structure in the resonance region 
gradually made progress [25MT] , the Regge analysis of the 
TTTT scattering amplitude was for a long time practically 
abandoned. 

The Regge parametrization used in [2H1] is based on 
Pennington's analysis [21]. This input was subject to a 
criticism in [32) , where an alternative Regge parametriza- 
tion was proposed. As shown in [33] , the driving terms and 
the 5* and P-phase shifts below 0.8 GeV are practically in- 
sensitive to the difference between the Regge parametriza- 
tion used in [5] and that proposed in [22- However, for 
extending the Roy analysis up to 1.15 GeV, an update of 
the Regge parametrization used as input is necessary. Ir- 



^ The maximum energy can be increased by using, instead of 
fixed-t dispersion relations as in the standard approach, more 
general dispersion relations along curves in the Mandelstam 
plane iSE]. 
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respective of the very modest impact it has on the region 
below 0.8 GeV, credit is due to Pelaez and Yndurain [32] 
for pointing out the need of improving the Regge repre- 
sentation used in [2]. In the meantime, these authors and 
their collaborators gradually improved their own Regge 
parametrization |34H37| . In the following we will compare 
our results with the most recent version of this series, [37], 
more precisely with the variant referred to as "CFD", 
which appears to represent the net result of their work. 

After a brief introduction to the Regge representation 
in the next section, we review the phenomenology of the 
leading Regge trajectories and of the total NN and ttN 
cross sections and discuss the consequences of factoriza- 
tion for the total tttt cross sections. The Regge represen- 
tation used in the present work is specified in section HI 
Section [5] contains a short description of the partial wave 
representation used in the present analysis, and section 
[S] presents the resulting picture for the total cross sec- 
tions. In section [71 we introduce the quantitative formu- 
lation of duality used in the present work and section |5] 
discusses the sum rules we are relying on to calculate the 
^-dependence of the Regge residues. Details concerning the 
way in which our analysis is carried out are explained in 
section in Our results for the Regge residues at t = are 
presented in section [TUl while those concerning their t- 
dependence are discussed in section [TTJ A few consistency 
checks are described in section [121 Finally, in section [T3l 
we present a summary and draw some conclusions. 

The paper has four appendices: in appendix[Xlwe give 
a short review of the Regge parametrizations of the tttt 
amplitude proposed in the literature. In appendix [Bl we 
specify the input used when solving the Roy equations for 
the partial waves (a detailed presentation will be given 
in a forthcoming paper [38] ) . In appendix [O we provide 
the explicit expressions of the sum rules used in this anal- 
ysis and in appendix [Dl we specify some details of the 
minimization procedure used to impose our constraints. 
Preliminary results of this work were reported in [39H43j . 



2 Regge poles 

The high energy behaviour of soft hadronic and photon- 
induced reactions can be described in terms of Regge poles 
and cuts in the angular- momentum plane. A Regge pole 
actually corresponds to the exchange of an infinite series 
of particles or resonances with quantum numbers such 
that they can be produced in the i-channel. For a de- 
tailed account of this framework and its application to 
meson-nucleon and nucleon-nucleon scattering we refer to 
the older literature |44f[51] . More recent studies [55H75] 
consider the description of the Pomeron in QCD [SB] and 
applications of Regge theory to difFractive phenomena and 
deep inelastic ep scattering at low x (see the excellent 
review [5T]). For soft hadronic processes, global fits of 
the ttN, KN and NN high energy data sets based on 
the Regge representation were performed in a system- 
atic way [62] , as quoted in the Review of Particle Physics 
[B31IMI- More sophisticated representations are obtained 
by including Regge cuts to describe absorption [651166) . or 



by the eikonalization of Regge exchanges in the impact- 
parameter representations [67H70] . These models extend 
the validity of the framework to larger momentum trans- 
fers and to higher energies, respectively. 

The contribution of a Regge pole to a scattering ampli- 
tude at large cm energy squared s and small momentum 
transfer t < Q has the form 



T{s,t) = -(i(t) 



exp[— i7rQ!(i)] - 



sm 7ra 



T I s 



'(«) 



(1) 



where a{t) and denote the trajectory and the residue, 
respectively, and r is the signature, taking the value 1 (—1) 
for C-even (odd) trajectories. The scale factor Si is usually 
set equal to 1 GeV^ and we will adopt this convention. Like 
the scattering amplitude T{s,t) itself, the residue /3(i) of 
a Regge pole is a dimensionless quantity. 

The trajectories a{t) and the residues /3(i) are ana- 
lytic functions of t, with branch points at the unitarity 
thresholds in the i-channel. The trajectories are univer- 
sal, i.e. a definite Regge pole in the angular momentum 
plane is located at the same place J — a{t) in all processes 
where it contributes. Information on the trajectories is in- 
ferred from the plot of the spin of known resonances ver- 
sus the square of their mass (Chew-Frautschi diagrams). 
Throughout this paper, we assume that in the t region of 
interest the trajectories are approximately lineai|3 



a{t) = a{0) + a'{0)t. 



(2) 



The residues /3(i) are much less known. The zeros of 
the denominator in (|T]) give rise to poles of the amplitude 
- unless the numerator or the residue /3(i) vanish there. 
For positive t, the poles correspond to the physical res- 
onances situated on the Chew-Frautschi plots. Poles in 
the spacelike region are unphysical, hence the residue /3(i) 
must vanish if a(t) -I- ^ (1-1- r) is an odd integer and t < 0. 
In the literature, this condition is sometimes implemented 
with factors depending explicitly on a(t) (see appendix 

A remarkable property of the residues, which follows 
from unitarity, is their factorization as a product of two 
vertices. For the exchange of a Reggeon R in the process 
a + b ~^ c + d, factorization reads 



f^nit) = lRac{t)lB,bd{t)- 



(3) 



For particles with spin, this relation holds separately for 
amplitudes of definite helicity. 

In the representation quoted by the Particle Data Group 
[63| , the dominant contribution, the Pomeron, is described 
as a triple pole in the angular momentum plane, which 
leads to an increase with energy as the square of a loga- 
rithm. Also, at nonleading order, daughter poles as well 
as Regge cuts arising from the simultaneous exchange of 
two or more poles contribute [49l[53]. In fact, in exotic 
channels, where resonances do not occur, the dominating 
contributions stem from such cuts. 



Deviations from linearity were also investigated, see for in- 
stance 1591. 
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The position acut{t) of the first branch-point arising 
from the exchange of two Regge poles is determined by the 
corresponding trajectories. For the exchange of two iden- 
tical poles, characterized by the linear trajectory a(0) -I- 
a'(0) t, the expression reads [i^[55] 

acut{t)=2a{0)-l + ^a'{0)t. (4) 

At large s, a Regge cut contributes as s"<^"*(*), with addi- 
tional logarithmic factors like ln'''(s/sc), where the scale 
Sc and the exponent 7 of the correction are in general 
unknown [^. Factorization does not hold for cuts and 
generally fails for absorption corrections |52| . especially 
at large negative t. 

3 Phenomenology of leading trajectories and 
total cross sections 

3.1 TvN and NN total cross sections 

During the last decade, a comprehensive compilation of 
the hadronic total cross sections was developed by the 
Compete Collaboration The Regge fit of the data 
above 5 GeV, quoted in PDG 2011 [53], is based on the 
expressions: 

<Jab{s) = Bln^is/so) + Zab+ (5) 
aab{s) = Bln^(s/so) + Zab + 

Yiabisi/sy^ +Y2ab (si/s^, 

where the target particle b is either a proton (p) or a neu- 
tron (n), and the projectiles a (a) are p or tt^, or their 
antiparticles. 

The first two terms on the right hand side of ([5]) rep- 
resent the contribution of the Pomeron. The coefficient B 
and the scale sq of the logarithmic term are assumed to 
be universal: 

B = 0.308 ±0.010 mb, = 5.38 ± 0.50 GeV . (6) 

The universality hypothesis was tested independently in 

izain]. 

The coefficients accounting for the energy independent 
part of the Pomeron contribution to the cross section are 

Z„„ = 35.63 ±0.25 mb, Z^+p = 20.86 ± 0.03 mb, (7) 

where the first value is the mean of Zpp and Zp„ (which 
agree within errors, indicating that the isospin breaking 
effects are very small). The numbers (H])-© show that - 
unless the energy is taken so high that ln(s/so) becomes 
large compared to 1 - the Pomeron term is approximately 
independent of the energy. 

The remaining terms in ([5]) account for the leading 
Regge poles, which for NN scattering are known to be 
/, 02, p and w [S3]. The connection with ([T]) is provided 



by the optical theorem, which relates the total cross sec- 
tions to the imaginary parts of the corresponding forward 
amplitudes. At high energies, it yields 

.po,e(.)«^(-j . (8) 

For NN scattering, the term proportional to Yi in ([5]) 
collects the contributions of the dominant C-even trajec- 
tories, / and 02, with the intercepts taken to be the same, 
'^/(O) — CKa2(0) = ^ ~ Vij the contributions from lower 
trajectories being neglected. Likewise, the leading C-odd 
trajectories, p and uj, are represented by the term propor- 
tional to Y2, with ap(0) — au,{0) — 1 — ?72- For complete- 
ness we quote the values from |63j : 

Yipp = 42.53 ±0.23 mb. Fa pp = 33.34 ± 0.33 nib (9) 
Yip„ = 40.15 ± 1.59mb, Fipn = 30.00 ± 0.96 mb. 

Using ([S|) for the collisions of protons or antiprotons 
on protons and neutrons, we can separate the contribu- 
tions from the individual trajectories, /, 02, p and w: 
Yipp — Yfj^jyj ± Ya^NN and Y2pp = 5^ivw ± YLfww- If thG 
proton target is replaced by a neutron, the isospin odd 
contributions from p and 02 change sign. Hence, the com- 
binations relevant for the scattering on neutrons are given 
by Y\pji — Yfj^j^ Ya2NN and Y2pn = ^w^v ± ^^w^v- 

At high energies, the cross sections barely show any 
difference between proton and neutron targets, indicat- 
ing that Ya^NN and YpjvN are small. In early works [46], 
the p and 02 contributions to pp and pp scattering were 
found to be negligible and in some analyses they were sim- 
ply dropped [37] • The recent Regge analysis of nucleon- 
nucleon scattering above 3 GeV reported in [751 confirms 
that / and uj dominate over 02 and p, respectively. Indeed, 
the values obtained from 

= 41.3 ±0.8 mb, Faajviv = l-2±0.8mb, (10) 
Yp^^ = 1.7 ± 0.5 mb , = 31.7 ± 0.5 mb 

show that the couplings of 02 and p to the nucleons are 
too small to clearly stick out. 

If isospin breaking effects are neglected, only the / 
and p trajectories are present in ttN scattering, in accord 
with the fact that two pions in a configuration with even 
(odd) isospin must carry even (odd) angular momentum, 
so cannot couple to 02 or oj. The parameters Y given in |63j 
for this process are: 

= 19.24 ±0.18 mb, Fp^^ = 6.03 ± 0.09 mb. (11) 
3.2 Leading trajectories 

The intercepts of the leading trajectories can be read off 
from the powers of s occurring in the representation ([S]) 
for the total cross sections. As mentioned above, this rep- 
resentation describes the Pomeron as a triple pole with 

ap(0)-l. (12) 
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Using the values r]i = 0.458±0.017 and 772 = 0.545±0.007 
from |63j and enlarging the errors to take into account 
other results in the literature (see appendix we take: 

a/(0) = 0.54 ±0.05, ap(0) = 0.45 ± 0.02 . (13) 

As stated above, we assume that in the t region of 
interest, the trajectories are adequately approximated by 
the linear formula ([2]). For the slope of the Pomeron tra- 
jectory, we use the value 

ap(0) =0.25 ±0.05 GeV~^ (14) 

which is consistent with the information on the shape of 
diffraction peak at high energies (see appendix El) . 

For af(t) and ap(t), we fix the slopes with the re- 
quirement that the trajectories pass through the masses 
squared: af{M^^) = 2 and ap{Mp) = 1, respectively. Us- 
ing the experimental values for the masses and adding 
errors to cover other results quoted in appendix [XJ we 
take 

a'^(O) = 0.90 ±0.05 GeV-^ (15) 
a'piO) = 0.91 ± 0.02 GeV"^ 

3.3 Factorization 

We rely on the values derived in section 13.11 in order to 
extract the Regge parameters of the mr amplitude through 
factorization, using equations ^ and ([5]) . As remarked in 
section [21 if the particles carry spin, factorization applies 
to amplitudes of definite helicity. In the forward direction, 
only the non-flip coupling survives, and factorization holds 
for each Regge pole contributing to the tttt, TriV and NN 
cross sections [53] . 

In the case of tttt scattering the representation ^ 
reads 

cr^±^+(s) = Sln^(s/so) + Z^TT + (16) 
Yfrr. (si/s)"^TV- (si/s)"^- 

The expression involves the leading Regge terms P, f and 
p. For the Pomeron, factorization gives the universality of 
the parameters B and sq and 

Z^^ = = 12.2 ±0.1 mb. (17) 

NN 

This is to be compared with the number obtained from 
"Z-KN with the quark counting rule: Z^tt ~ | ■^ttn ~ 14 mb. 
For the Y parameters, the values in (jl0|) and dTl]) give 

yl 

y^^^ = =8.95 ±0.24 mb, (18) 

^ f NN 

Yp^^ = = 21.8 ± 9.0 mb. 

^ pN N 

The phenomenological information about the coefhcients 
Z and Y will be converted into standard Regge terminol- 
ogy in section 



3.4 Data on tttt cross sections 

We complete this phenomenological discussion with a brief 
review of the experimental results for the tttt total cross 
sections obtained indirectly, from reactions like ■K~p — >■ 
TT+TT-n, TT-p TT--K-A++ ^MW\ and TT^p 
TT^n^pX [12113]. 

Fig.[T]shows that in the 7r~7r+ channel, the three avail- 
able data sets [7i .9al2) are more or less consistent, not only 
with one another, but also with the cross section that fol- 
lows by applying factorization to the Regge representation 
in PDG 2011. The figure also indicates the cross section 
corresponding to some of the Regge representations avail- 
able in the literature, as well as the outcome of our own 
analysis. These entries will be discussed in section [10.51 




Vs (GeV) 

Fig. 1. Total TT^-K^ cross section 

In the channel tt^tt*^, we are aware of only one set of 
data [7]. As these concern the region below 2 GeV, they 
do not contain information about the behaviour at high 
energies. 

The TT^vr^ channel does not contain resonances. At low 
energies, the cross section is small and notoriously difficult 
to measure. Various data are shown in Fig.[51 Those of 
Cohen et al. [10] concern the elastic rather than the total 
cross section. At low energies, however, the inelasticity in 
this channel is very small, so that the elastic cross section 
barely differs from the total one. In particular, the main 
source of inelasticity in 7r~7r+ collisions - transitions into 
KK states - is forbidden here. This implies that the data 
of Biswas et al. [7] and Losty et al. [11, are in conflict 
with those of Cohen et al. [TU] . Concerning the behaviour 
at higher energies, there is a clash between the data of 
Robertson et al. [S] and Abramowicz et al. [13] . Note that 
the energy range of this plot extends to 20 GeV, far beyond 
the region which matters in our context. 

In the Regge parametrization of the tt^tt^ cross sec- 
tion, the p contributes with a sign opposite to Pomeron 
and /. Since factorization determines the latter rather 
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Fig. 2. Total n n cross section 



sharply, the uncertainty of the result for cr^-j^- is dom- 
inated by the one in the residue of the Regge pole as- 
sociated with the p. The red shaded area indicates the 
result obtained from ttN and NN scattering with factor- 
ization, using equations (jl7p and (llSp . The uncertainty in 
this result is large, because (a) the p contribution to NN 
scattering is small and (b) the extraction of the p residue 
requires data on p n scattering, which are less precise. The 
other entries shown in the figure will be discussed in sec- 
tion [1031 



4 Regge representation of the tttt amplitude 

We now specify the parametrization of the tttt scattering 
amplitude we are going to use at high energies. The dom- 
inant Regge poles with C — 1 are the Pomeron and the 
/, while the channel with C = — 1 is dominated by the p- 
pole. What counts is the isospin in the i-channel. We use 
the notatior0 of reference [2] , where the component of the 
scattering amplitude with i-channel isospin It is denoted 
by T^^*\s,t). In this notation, the optical theorem takes 
the form 

V-'^(s-4M2) 

The Pomeron and the / contribute to T^^"^ and p con- 
tributes to T^^\ while the exotic amplitude T^^^ only re- 
ceives subdominant contributions from the exchange of 
two or more Regge poles. In the dispersive framework we 
are using, the real parts are represented as dispersion in- 
tegrals over the imaginary parts. As we are evaluating 
these integrals only at low energies, it suffices to specify 
the Regge representation of the imaginary parts - for the 
analysis of the Roy equations and sum rules discussed in 



^ In [321137] the scattering amplitude is instead denoted by 
F{s,t). The normalization differs from ours by a numerical 
factor: T{s,t) = 4tv^ F{s,t). 



the present paper, the high energy behaviour of the real 
parts is not relevant. 

The quantities cr^±^+ , <J^±„a (s) and u^o^o (s) represent 
linear combinations of the total cross sections of definite 
isospin in the t channel: 

a.±,+ (,s) ^ i (,s) T I (s) + i ^(2) (,) ^ (20) 



f^7r±7ro(s) 



r(") 



2 

1 ^(2) 
3 



2^(2) 
3 



(.). 



4.1 Parametrization of IniT^°\s,t) 

Our explicit representation for the imaginary part of the 
component with It = reads: 



ImT(0)(s,t)R, = /3p(t) - 



ap{t) 



(21) 



1 + B In^ ( - 

So 



af(t) 



It contains the trajectories of Pomeron and /, 

ap{t) = ap{0)+a'p{0)t, 
af{t) = af{0)+a'f{0)t. 



(22) 



The range used for the intercepts and slopes is specified 
in equations (IT^ - p^ . 

We first comment on the coefficients relevant at i = 0, 
which govern the total cross section (T^°)(s). Using the 
optical theorem (|19p on the right hand side of the first 
equation in ([20)1 and comparing the result with ([T5|. we 
obtain: 



^p(O) = 3 siZ^^ihc)-^ = 94 ± 1 , 
B = B/Z^^ = 0.025 ± 0.001, 



(23) 



where we have explicitly indicated the conversion factor 
{hey = 0.389 mbGeV^ needed to arrive at dimensionless 
residues. As discussed in appendix 1X1 the values of /3p(0) 
and /3/(0) proposed in the literature cover a large range. 
The relative weight of the contributions from the Pomeron 
and from the / depends on the parametrization chosen for 
the Pomeron. For our analysis, this does not present a sig- 
nificant source of uncertainty: the relative weight of the 
two contributions barely matters. It is important, how- 
ever, that in the energy range of interest, factorization 
yields an unambiguous prediction for the sum of the two, 
that is for the total cross section with It — 0. 

In (HI]), we have allowed for a pre-asymptotic contri- 
bution, whose size is determined by the coefficient pi . The 
reason for including such a term is that we are making use 
of the parametrization in terms of Regge poles down to an 
energy of 1.7 GeV. The uncertainties in the residues given 
in (|23l) stem from an analysis of the observed ttN and 
TV TV cross sections above 5 GeV. They do not cover the 
contributions generated by nonleading Regge poles, which 
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become increasingly important as the energy is lowered. 
In the representation (|2T|) . these contributions are mod- 
eled by a single term, which drops off more rapidly with 
the energy than the Pomeron, by one power of s. 

We now turn to the ^-dependence. While the term as- 
sociated with the / is of the standard form of a Regge pole 
given in equation ([1]), the one describing the Pomeron con- 
tains the square of a logarithm. Following |63j . we made 
the assumption that the Pomeron is a triple pole in the 
angular momentum plane, located at J = ap{t). This is 
legitimate only if the energy is not too high. As pointed 
out in [ZnilTj, unitarity imposes a bound on the shape 
of the diffraction peak. The parametrization in (j2ip vio- 
lates this bound at very high energies, for the following 
reason. In our normalization of the scattering amplitude, 
the differential cross section for elastic scattering reads 



da _ \T(s,t)\^ 

H ~ 167rs(s-4M2) ■ 



The contribution to this cross section from |ImT(s,t)p 
alone grows in proportion to In^s s^"^''-"-*', so that the in- 
tegral over t - the total elastic cross section - increases 
with the third power of In s and thus eventually becomes 
larger than the total cross section. The same problem 
also occurs if the intercept of the Pomeron is taken above 
unity, ap{0) > 1. There are well-known ways to avoid the 
clash with unitarity, for instance in Regge-eikonal mod- 
els [57H7D] valid at very high energies. 

In the context of the present paper, the behaviour at 
very high energies is an academic issue, because the inte- 
grals over the Regge representation occurring in the eval- 
uation of the various quantities of interest pick up signif- 
icant contributions only from modest energies, where the 
logarithmic term represents a small, roughly energy in- 
dependent correction. It could just as well be absorbed 
in the residue Pp{t) and in the sub-leading coefficient 
Pi, so that the problem would disappear. Nevertheless, 
a parametrization of the Pomeron that (i) reduces to the 
one used in the Review of Particle Properties [63] when 
t — ^ and (ii) is consistent with general principles also for 
large s and t would be of considerable interest - we 
did not find such a representation in the literature. 

We give the term proportional to pi the same t-depen- 
dence as the Pomeron. This has the advantage that the 
pre-asymptotic contributions in the channel with It = 
are accounted for in terms of a single parameter, but it is 
evident that a simple parametrization of this sort can ac- 
count for the nonleading terms only in a very crude way. 
Nevertheless, we will show that this parametrization al- 
lows us to make contact between the region where the first 
few partial waves represent a reliable approximation and 
the one where the amplitude can be represented in terms 
of a few Regge poles. The t-dependence of the residues 
I3p{t) and f3f{t) is discussed in section [TT] 



4.2 Parametrization of ImT'^^\s,t) 

For the component with = 1, we work with the repre- 
sentation 

Imr(i)(s,t)R = /3p(t){l + ri^} (^^^ . (25) 

It involves the p-trajectory, 

ap{t) = ap{0) + a'p{0)t, 

whose parameters are specified in (jl3l) and (jTSj) . as well 
as a pre-asymptotic contribution, determined by the pa- 
rameter ri. For the intepretation of this term, we refer 
the reader to the last paragraph of the previous section - 
what we wrote for pi also applies to ri . 

The result for the coefficient YpT^^^ in ((HI), which fol- 
lows from the Regge fits to the ttN and NN data in |63) . 
comes with a large uncertainty. The same applies to the 
corresponding result for the p residue at t = 0: 

/3p(0) = 2siV,(?ic)-2 = ii2± 46. (26) 

In the following, we do not make use of this crude estimate. 
Instead, we rely on the determination of (3p{0) based on 
the Olsson sum rule [IS] , which will be discussed in section 
110.21 Our determination is consistent with the above num- 
ber, but somewhat more precise. The behaviour of /3p(i) 
at nonzero values of t is discussed in section 1111 

4.3 Parametrization of Imr(^)(s,*) 

At the energies considered in the analysis of the high en- 
ergy ttN and NN data, contributions from exotic Regge 
channels are neglected. In the energy range discussed in 
the present paper, however, this is not justified. Indeed, as 
we will see, some of the sum rules obeyed by the scatter- 
ing amplitude do require a small component with It — 2. 
In this channel, the dominating contributions stem from 
the exchange of two Regge poles with /t = 1, so that the 
leading singularity in the angular momentum plane is a 
branch cut |21j . The leading contributions stem from p—p 
and 02 — 02 exchange. In the Regge representation we are 
using, which is based on [53] , the 02 trajectory is taken de- 
generate with the /. Moreover, equations (jl3L (|T5|) show 
that the / and p trajectories are rather close to one an- 
other: for the cut generated by p—p exchange, the formula 
gives app(O) = - 0.10 ± 0.04 and a;,p(0) = 0.46 ± 0.01, 
while the cut from / — / exchange is characterized by 
a//(0) = -h0.08±0.10, a^^(O) = 0.50 ± 0.08. Also, at 
the relatively low energies considered here, the neglected 
logarithms cannot be distinguished from the difference be- 
tween app(O) and a//(0). The numerical values indicate 
that, compared to the dominating contribution from the 
Pomeron, the exotic component of the amplitude roughly 
falls off in proportion to 1/s. The It = 2 component of 
the scattering amplitude may thus be viewed like a pre- 
asymptotic contribution. We approximate it in the same 
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manner as the one occurring in the amplitude with It — 0, 
parametrizing it with a simple Regge pole, 

Imr(2)(s,i)R = /3eW(^^) (27) 

where (3e(t) represents an effective residue and 

ae(t) = ae(0) + a^(0)i (28) 

is an effective trajectory with ae(0) — 0. 

The formula ^ implies that the slope of the effective 
exotic trajectory is smaller than the slopes of the p or /- 
trajectories, by a factor of 2. This indicates that, compared 
to p and /, the relative importance of the exotic amplitude 
increases as one goes away from t = 0, towards spacelike 
values [61]. On the other hand, according to (fT4)) . the slope 
of the Pomeron trajectory is even smaller, so that, com- 
pared to the Pomeron, the exotic component loses weight 
if t becomes negative. We fix the exotic trajectory with 

ae(0)=0, a'g(O) = 0.5±0.1GeV~^ (29) 

The properties of the residue l3e(t) as a function of t are 
discussed in section [TT] 

5 Partial waves 

The representation of the scattering amplitude in terms of 
Regge poles and cuts is dual to the one in terms of partial 
waves: the amplitude can be represented with either one 
of the two. We denote the partial wave amplitudes by 
t({s), where / = 0,1,2 stands for the s-channel isospin 
and ^ = 0, 1, 2 . . . denotes the angular momentum. We use 
the normalization 

p{s) = v/l-4M2/s, 

where 77^ (s) and Sj{s) denote the elasticity and the phase 
shift, respectively. The partial wave decomposition of the 
scattering amplitude then reads 

T'is,t,) = 32^^(2^+l)P,(z)4(s) (31) 

e 

= i(4M2-s)(l-z), 

where Pe.{z) is the Legendre polynomial of order For the 
even angular momenta, Bose statistics permits / = and 
/ = 2, while waves of odd angular momentum carry / = 1. 
We use the symbols S° and to distinguish the two S- 
waves, write P^ for the P-wave and label the higher partial 
waves analogously. In the region where the partial waves 
are small, the real and imaginary parts are approximately 
related by 

luitlis) ^ p{s)Ket\{sf + . (32) 

2 pis) 



In this region, unitarity (0 < ?/^(s) < 1) sets a lower bound 
to the imaginary part, given by the first term in the above 
relation. 

At low energies, the angular momentum barrier sup- 
presses the contributions from high angular momenta. The- 
re, the decomposition of the amplitude into partial waves 
is preferable, because the first few terms in this decompo- 
sition suffice to obtain a decent approximation. At high 
energies, on the other hand, the Regge representation is 
more efficient, because the high energy behaviour of the 
amplitudes is dominated by the leading Regge poles. In 
the context of the present paper, which aims at an im- 
proved Regge parametrization, the partial wave represen- 
tation plays an important role for two reasons: (i) We are 
making use of sum rules to analyze the ^-dependence of 
the Regge residues. Since these involve integrals extending 
from threshold to infinity, we need an adequate represen- 
tation also at low energies, (ii) The leading Regge poles 
yield a good approximation for the scattering amplitude 
only at very high energies. We need to bridge the gap be- 
tween the region where the leading Regge poles dominate 
and the one where the Roy equations provide a reliable 
representation for the partial waves. Since the behaviour 
of the amplitudes between these two domains represents 
the main source of uncertainty in our analysis, we need to 
discuss it in some detail. 

In the region where the Roy equations are valid, we 
determine the partial waves by solving these equations. 
The input of the calculation consists of four parts: 

1. subtraction constants 

2. elasticities below Smax 

3. partial wave imaginary parts above Smax 

4. driving terms 

where ^/Smax is the upper end of the energy interval over 
which the equations are solved. In [JJ,*?!, the upper end 
was taken at ^Smax =0.8 GeV and the analysis was re- 
stricted to the S- and P- waves. In that case, the input 
unambiguously determines the solution. In the present pa- 
per, we consider the entire range where the Roy equations 
are valid, i.e. set ^Sraax — \/68M^ ~ 1.15 GeV, and also 
calculate the D-, F- and G- waves. As discussed in detail 
in [2] , the input listed above does then not determine the 
solution uniquely: on the extended interval, the Roy equa- 
tions admit a 3-parameter family of solutions. We find it 
convenient to identify two of these degrees of freedom with 
the value of the S° phase shift at two points in the inte- 
rior of the interval on which we solve the Roy equations: 
5^(sa), 8l{^M\), with = 0.8 GeV. The third one is 
fixed with the P^ phase shift: 8\[sa). The phenomenolog- 
ical estimates used for these quantities as well as those for 
the elasticities and for the imaginary parts of the partial 
waves are discussed in appendix IB. 31 The calculation of 
the driving terms is sketched in section [HI 

Near threshold, the subtraction constants control the 
behaviour of the scattering amplitude. In the case of tttt 
scattering, the dispersive representation involves two sub- 
traction constants, which may be expressed in terms of 
the two S-wave scattering lengths, Og, Og. At leading or- 
der of the chiral expansion, their values can be predicted 
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in terms of the pion decay constant [75]. The higher or- 
der corrections have been calculated in the framework of 
xPT [TliiSOj. As pointed out in [3], the dispersive analy- 
sis of the scalar pion form factor |5T] can be used to pin 
down one of the low energy constants occurring therein, 
so that a remarkably sharp theoretical prediction for the 
scattering lengths follows: 



0.220 ±0.005, 



2 
Co 



-0.0444 ±0.0010. 



(33) 



In the meantime, this prediction has been tested in a series 
of beautiful low energy precision experiment^ concerning 
the decays K Hv-ktt [83i^, K — >• tttttt [85] and pio- 
nic atoms [55]. For the combination aj] — Og, where the 
experimental precision is comparable to the accuracy of 
the prediction, theory and experiment agree within errors. 
The data do not pin down the two individual scattering 
lengths to the same accuracy, but in combination with the 
fact that the dispersive evaluation of the scalar form factor 
leads to a sharp correlation between Aq and Oq, the exper- 
imental results imply = 0.2198(46)stat(16)syst(64)t,j, 
al = -0.0445(ll)stat(4)s;/st(8)t/i, thus confirming the pre- 
dictions to a remarkable degree of accuracy. Alternatively, 
as pointed out by Stern and collaborators (HJ, the two 
scattering lengths may be disentangled by invoking com- 
paratively crude experimental results extracted from TriV 
scattering. The recent update of this analysis by Garcia- 
Martin et al. [771 leads to ajj = 0.220(8), al = -0.042(4), 
also consistent with the predictions, albeit somewhat less 
precise. 

Moreover, the recent progress made on the lattice now 
allows a calculation of Oq from first principles, analyzing 
the volume dependence of the energy levels of the system 
in a box of finite size [55^150] . Also, the low energy con- 
stants ?3, ?4, which dominate the uncertainties in the the- 
oretical prediction for aj] and Oq, can now be determined 
on the lattice, from the dependence of M.^ and F.^ on the 
quark masses. The results are consistent with the esti- 
mates used in [3] and the uncertainties in 1^ are reduced 
significantly [HI]. All of this corroborates the conclusion 
that the subtraction constants are reliably and accurately 
known. The remaining three parts of the input will be 
discussed below. 



6 Partial wave representation of the total 
cross sections 

In the present section, we briefly describe the qualitative 
picture obtained for the total cross sections, in order to 
set the stage for our analysis. The method used to derive 
this picture is outlined in section [H| 

We denote the total cross sections with s-channel iso- 
spin Is — 0,1,2 by cr'^(s), cr^(s), cr^(s), respectively. In 
this notation, the total cross sections for pions of definite 
charge are given by 



C7,+ ,+ (s) =CT^-,-(s) =ct'(s) 



(34) 
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the total cross sections obtained from the sum over the 
partial waves. The dashed black lines represent the central 
solutions of the Roy equations for the S- and P-waves. The 
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contributions from these waves dominate the total cross 
sections below 1 GeV, reach a maximum and then drop 
with the energy. The dotted black curves are obtained by 
adding the imaginary parts of the central Roy solution 
for D° and D^, respectively, while the dash-dotted ones 
include the F- and G- waves. 

In the total cross section, the resonance peaks occur- 
ring in the imaginary parts of the individual partial waves 
generate an oscillatory behaviour, which only disappears 
gradually with increasing energy. Fig. 3a shows that, below 
1 GeV, cr'^(s) is dominated by the contribution from the 
S^-wave, which exhibits a broad bump around 600 MeV. 
The phenomenon reflects the presence of a pole on the 
second sheet of the partial wave amplitude to(s), the well- 
known /o(600) or a. Slightly above 0.8 GeV, the phase 
shift passes through 90°. There, the inelasticity is very 
small, so that the S°-wave nearly reaches the unitarity 
limit. This explains why the uncertainties in the total cross 
section are very small in that region. While the destructive 
interference with the /o(980) generates a pronounced dip 
in the vicinity of KK threshold, the behaviour between 
1.15 and 1.5 GeV is dominated by the /2(1275). The con- 
tribution from the G^-wave becomes visible only above 
1.6 GeV. The first resonance with these quantum num- 
bers, the /4(2018), can be seen, but it does not produce a 
marked peak, because many decay channels are open and 
the branching fraction r{f4 7r7r)/P(/4 — >■ all) is below 
20%. 

The most prominent resonance, the p{770), manifests 
itself in (s). Since an appreciable inelasticity sets in only 
at + AU ^ 0.92 GeV, the contribution from the P^- 
wave nearly saturates the unitarity limit at the peak of 
this resonance. The cross section a^{s) exceeds 220 nib 
there - the upper two thirds of the peak are chopped off, 
to make the behaviour of the cross section above 1 GeV 
more visible. The spin 3 resonance p3(1690) produces a 
more modest peak, followed by an even weaker structure 
around 2 GeV. In the PDG listings, the lowest resonance 
with ^ = 5 is the p5(2330). We expect partial waves with 
£ > A to play a significant role only above 2 GeV. 

The exotic cross section cr^(s) — cr^-^- (s), which does 
not receive contributions from resonances, is shown in 
Fig. 3c, a magnified version of Fig. 2. Below 1.15 GeV, the 
Roy equations provide good information about this cross 
section: they imply that the real parts of the exotic par- 
tial waves are dominated by the contributions from the 
imaginary parts of the non-exotic waves, which are known 
comparatively well. The uncertainties in cr^(s), however, 
grow with the energy - at ^/s = 2 GeV, they are of order 
3 mb. 

Fig. 3 shows that, in the region between 1.5 and 2 GeV, 
our partial wave representations for the total cross sections 
with Is — and Is — 2 are roughly energy independent, 
while the one with Is = 1 clearly reveals the presence of 
a resonance. Our quantitative estimates for the imaginary 
parts of the partial waves, which are outlined in appendix 
IB.3[ lead to substantial uncertainties in the total cross sec- 
tions, which, moreover, grow with the energy - beyond 2 
GeV, the partial wave decomposition is barely of practi- 



cal use in our context. In the channels with = or 1, 
the uncertainties mainly stem from the background under- 
neath the resonances. In this language, the cross section 
in the channel with Is = 2 represents pure background. 
Indeed, Fig. 3 shows that our estimate for the uncertain- 
ties in cr*'(s) and a^{s) are comparable to the exotic total 
cross section cr^(s). 



7 Duality 

In the preceding sections, we discussed the Regge and par- 
tial wave representations by themselves. We now merge 
the two and first consider the forward direction, t = 0. 

A simple way to join the two parametrizations of the 
imaginary parts is to use the partial wave representation 
below some energy Et and the Regge representation above 
that energy. In order for the total cross sections not to 
make a jump at ^/s = Et, the imaginary parts must be 
continuous there: 

• In the case of ImT'^''^(s, 0), this requirement can be 
used to fix the coefficient pi of the pre-asymptotic term, 
so that, at i = 0, the Regge representation (|21l) does then 
not contain any unknowns: the total cross section cr'^'^^ (s) 
is determined also above the transition energy Et- 

• The Olsson sum rule subjects the cross section ct^^^ (s) 
to a constraint, which correlates the value of /3p(0) with 
the pre-asymptotic coefficient ri (see section [10.2p . Since 
the continuity requirement yields a second constraint of 
this type, we can determine /3p(0) as well as ri. As a re- 
sult, the total cross section is determined above Et , also 
in the channel with /( = 1. 

• Finally, in the exotic ^-channel, continuity deter- 
mines the residue /3e(0) in terms of the partial wave repre- 
sentation at Et, so that the Regge representation of cr^^-' (s) 
does then not contain any free parameters, either. 

The choice of the transition energy represents a com- 
promise: if Et is taken too low, then our Regge represen- 
tation does not provide a decent approximation, if it is 
taken too high, then our neglect of the partial waves with 
^ > 4 is not justified. 

In [2]-[4] , the transition from the partial wave represen- 
tation to the Regge representation is made at 2 GeV, while 
in |5^[511 - I57] the Regge parametrization is assumed to be 
valid down to 1.42 GeV. The problem with the former 
analysis is that the phenomenological information used 
for the partial waves above 1.4 GeV leaves much to be 
desired, while the problem with the latter choice is that 
Regge asymptotics does not set in that early - see Fig. 3b. 

Continuity may be viewed as a local version of duality: 
it requires the partial wave and Regge representations to 
agree at one particular energy. In the following, we use 
a weaker form of this requirement: we assume that the 
partial wave and Regge representations agree with one an- 
other only in the mean. For short, we refer to the resulting 
constraints as duality conditions. 

More precisely, we replace the continuity conditions 
discussed above by the following constraints. Fig. 4 shows 
that the total cross sections with fixed i-channel isospin 
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either exhibit a peak or a dip, generated by the p(1690). 
The dashed vertical hnes indicate the adjacent valleys at 
1.5 and 1.9 GeV, respectively. The bands labeled "Regge" 
are obtained by assuming that, on the interval spanned 
by these lines, 1.5 GcV < y/s < 1.9 GeV, the average 
over the partial wave representation of the total cross sec- 



tions agrees with the average over the Regge representa- 
tion. In other words, we assume that the inclusion of the 
pre-asymptotic terms allows us to extrapolate the Regge 
representation down to the above interval, not locally, but 
in the mean. Roughly, this amounts to imposing continu- 
ity at Et, varying the value of Et over the above range 
and averaging the result. 

This version of duality lends itself to a natural general- 
ization as well: it suffices to replace the total cross sections 
(T^'^^ (s) by the functions p{s) ImT^*'^ (s, t). The duality con- 
ditions used in the present paper represent the require- 
ment that the mean value of the Regge representation for 
these functions agrees with the partial wave representa- 
tion thereof. Duality is useful in our framework, because it 
provides us with an estimate for the pre-asymptotic terms 
Pi, ri, /3e(0), much like the continuity conditions discussed 
above. Moreover, its generalization provides a constraint 
on the i-dependence of the Regge residues. 

Whenever we calculate integrals over the imaginary 
parts of the tttt amplitude we need to specify how the 
transition between the partial- wave and the Regge repre- 
sentation takes place. We assume that the partial wave 
representation provides an adequate approximation below 
1.7 GeV and that the Regge representation yields an ad- 
equate approximation above 2 GeV. In the transition re- 
gion, we use an interpolation, 

ImT(^')(s,i) = {1- /i(s)}ImT(-f')(s,i)pw (36) 

+ /l(s)ImT(^*)(s,i)R, Sa<S<Sb, 

with — 1.7 GeV, = 2 GeV. The function h{s) in- 
terpolates between h{sa) = and h{st,) — 1. For definite- 
ness, we take a cubic polynomial and fix the coefficients 
with the requirement that the first derivative is contin- 
uous at the lower and upper ends of the interpolationlfl 
The interpolation p6p automatically ensures continuity, 
also at non-zero values of t. 



8 Sum rules 

The Olsson sum rule [TS] originates in the fact that the 
t-channcl 1=1 amplitude does not receive a Pomeron 
contribution, and thus grows only in proportion to s"p(*) 
for s — > CXI. The fixed-t dispersion relation obeyed by 
this amplitude, however, does contain terms that grow lin- 
early with s. For the relation to be consistent with Regge 
asymptotics, the contribution from the subtraction term 
must cancel the one from the dispersion integral. At i = 0, 
this condition reduces to a sum rule which expresses the 
combination 2ag — 5al of the S'-wave scattering lengths as 
an integral over the cross section (t^^^(s) |15) : 



2aO 



5ao = 



3Ml 
47r2 



4M2 vs(s 



4M^) 



(37) 



The requirement that the amplitude T^^\s,t) has the 
proper high energy behaviour also for t ^ implies a 

^ The corresponding explicit expression for the interpolating 
function reads h{s) = (s — Sa)'^(3 Sb ~ Sa — 2 s)/ (sj, — Sa)'^- 
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further constraint. This amphtude obeys the fixed-i dis- 
persion relation in equation 



66t of appendix IB.ll The 
constraint can be derived by evaluating the coefficient of 
the term on the r.h.s of this equation that grows linearly 
with s and requiring that the term has the same value as 
for t ^ 0. This condition subjects the function Im T*^^) (s, t) 
to an entire family of sum rules of the form 



(38) 



where t is a free parameter and the function S{t) is given 
by: 



(39) 



1 



ds 



AMI 



2(s-2M2)ImTi(s,0) 
s {s - 4M2) {s-t){s + t- 4M2) 



The first integral involves only the It — 1 component of 
the imaginary part. The bar indicates that the value at 
t = is subtracted: 



T{s,t) 



T{s,t)-T{s,0) 



(40) 



The function T^(s,0) in the second integral is the = 1 
component at i = 0, which exclusively picks up contribu- 
tions from partial waves of odd angular momentum. The 
integrand in (15^ displays a fictitious pole for s = AM^ — t. 
However, for a crossing symmetric integrand, the residue 
of the pole vanishes. The sum rule S{t) = was written 
down already in [TB] and was applied for the determination 
of the Regge residues in [T71IT5] . For a detailed discussion, 
we refer to appendix C in |2]. 

The Roy equations are manifestly symmetric with re- 
spect to s -H- M, but they do not ensure crossing symmetry 
with respect to s •(-> i, which requires the scattering am- 
plitude T = (T°,r\r2) to obey the relation 



Tis,t)-CstT{t,s)=0. 



(41) 



On the other hand, any solution of the Roy equations that 
obeys this constraint is then automatically also crossing 
symmetric with respect to t ^ u. Note that the subtrac- 
tion constants as well as the contributions from the S- and 
P-waves drop out here: the relation (j4ip amounts to a fam- 
ily of sum rules that relate integrals over the imaginary 
parts of the partial waves with £ > 2 to integrals involv- 
ing the Regge residues. We exploit these additional con- 
straints by evaluating the derivative of (|4T|) with respect 
to s at s = 0. As shown in appendix [C] this leads to a new 
family of sum rules of the same type as the one above. We 
write these in the form Ck{t) = 0, with fc = 0, 1, 2. Explicit 
expressions for the three functions Co{t), Ci{t), C2{t) are 
given in equation (j75l) . As we will see, these sum rules pro- 
vide us with a good tool to determine the ^-dependence of 
the residues. 



9 Method of analysis 

Our analysis runs as follows: 



1. We start with a complete representation of the ampli- 
tude, which specifies the imaginary parts of the partial 
waves below 2 GeV as well as the parameters occurring in 
the Regge representation. 

2. As a first step, we calculate the driving terms for the 
partial waves with £ < 4. The driving terms of a given 
wave consist of two parts: contribution from the imagi- 
nary parts of the other partial waves with £ < A below 2 
GeV and contribution from all waves above 1.7 GeV. The 
former is smoothly cut off with the factor {l — h{s)}, while 
the latter is gradually turned on with the factor h{s), in 
accord with equation ((36|) . 

3. Next, we specify the input for which we wish to de- 
termine the full amplitude. Some of the input variables 
concern the partial wave representation, others concern 
the Regge parametrization. In fact, neither one of the 
two representations is fully specified by the input vari- 
ables alone (see appendix IB. 3p . We refer to the remaining 
parameters as output variables. The Roy equations deter- 
mine the partial wave output, while the sum rules and 
duality conditions serve to determine the Regge output. 

4. Together with the driving terms, the partial wave input 
uniquely determines the solution of the Roy equations. 
We solve these equations iteratively. The result is a new 
complete representation of the amplitude. 

5. We then evaluate the various sum rules and duality 
conditions, which amount to constraints imposed on the 
Regge output variables. Since some of the sum rules and 
the duality conditions depend on t, we are in principle 
dealing with an infinite number of equations for a finite 
number of parameters and can therefore not expect to find 
an exact solution. Instead, we determine the minimum of 
a discrepancy function, obtained by adding up the various 
equations in quadrature ~ for details, we refer to appendix 
Id] The minimum fixes the output variables of the Regge 
parametrization. In combination with the input specified 
in step 3, we thus arrive at a new Regge representation. 

6. If the result of the preceding step differs significantly 
from the Regge representation we started with, we must 
return to step 2 and iterate the entire procedure. We find 
that the process converges rapidly - two or three itera- 
tions suffice to arrive at a self-consistent solution of good 
quality. 

7. For a given choice of the input, the above procedure pro- 
vides us with a complete representation of the amplitude. 
We finally need to determine the response of the output 
to changes in the input. As the variations of interest are 
small, the result for the partial wave and Regge represen- 
tations of the scattering amplitude is approximately lin- 
ear in these variations. We treat the input parameters as 
independent Gaussian variables, with one exception: we 
account for the correlatior0 between the two scattering 
lengths Oq and Oq. The error analysis is then straightfor- 



This correlation cannot be ignored; the error in the com- 
bination 2oo — Soq, for instance, which plays a prominent role 
in the Olsson sum rule, is significantly smaller than what is 
obtained if the two terms are treated as independent. 
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ward - the quantities of interest are linear functions of 
the input variables and we can calculate the uncertainties 
therein as well as the correlations between them by aver- 
aging over the Gaussian fluctuations in these variables. 



10 Regge analysis of the total cross sections 

10.1 Total cross section with It — 

We first illustrate the method with the determination of 
the pre- asymptotic coefficient pi. As discussed in section 
[21 factorization yields a rather sharp prediction for the 
asymptotic behaviour of the total mr cross section in the 
channel with It = 0. The prediction is shown as a narrow 
red band in Fig. 4a, obtained by inserting the values (^5]) 
for fip{0) and /3/(0) in the formula ((2T|) . setting the pre- 
asymptotic coefficient pi equal to zero and calculating the 
corresponding cross section with the optical theorem (|19p . 

Since the grey band in Fig. 4a, which represents the 
partial wave representation of the cross section a^^\s), 
passes below the extrapolation of the contributions from 
/3p(0) and /3/(0), the pre-asymptotic term must be neg- 
ative. The numerical evaluation of the duality condition 
formulated above implies that the coefficient pi is in the 
range 

Pi = -0.8±0.6. (42) 

10.2 Olsson sum rule 

The analysis of the total cross section in the channel with 
It — 1 is more involved, because we do not make use of 
the result obtained from the ttN, NN and NN data with 
factorization. As discussed in section [4.21 the error bar in 
that result is rather large. We instead rely on the Olsson 
sum rule ([57)) . which states that an integral over the cross 
section ct*^^^(s) is determined by a combination of S-wave 
scattering lengths. The l.h.s. is known from xPT [5]: 

2a[J-5a2 = 0.663 ±0.007. (43) 

The prediction is very sharp, because the NLO corrections 
to Weinberg's low energy theorem (78) exclusively involve 
the low energy constant €4 ^79] . Chiral symmetry implies 
that the same constant also governs the chiral expansion 
of the scalar radius of the pion. The value quoted in (|43l) 
relies on a dispersive calculation of the scalar pion form 
factor |81) and accounts for all contributions up to and 
including NNLO of the chiral expansion ^80j . As discussed 
in section [51 the predictions for the scattering lengths have 
in the meantime been corroborated, both experimentally 
and on the lattice. Concerning the specific combination 
relevant for the Olsson sum rule, the "best value" obtained 
on the basis of the available experimental information in 
123 is 2ag - 5al = 0.650 ± 0.015, thus confirming the 
prediction (|43p to an accuracy of 2.5%. 

The integral on the r.h.s. of (j37p can be evaluated with 
the optical theorem Below 1.7 GeV, the integrand is 
represented as a sum over the partial waves discussed in 



section [SI while above 2 GeV, the Regge representation in 
(|25l) is relevant. In the transition region, we are using the 
interpolation p6p. The Olsson integral is then given by 
the sum of two terms: O — Opw + Or, with 

op„.2if ,44) 

47r^ ^s{s - 4M2) 

_ 3M2 h{s) (1) 

Or = , „ / as — = ctr Is j . 

Below 1.15 GeV, we solve the Roy equations for the S, 
P, D, F and G partial waves. From there to 2 GeV, we use 
a parametrization based on phenomenology (see appendix 
IB.3P . Table [T] lists the numerical results obtained for the 
contributions from the partial waves considered. The table 



pi 

DO 
D^ 
pi 
GO 
G^ 


n o-|C,+0.019 

-U.U78_o 
246+0 01'^ 

U.Z'iO_o 019 

0.061 ± 0.003 
-0.007 ±0.003 
0.013 ±0.001 
0.0021 ± 0.0003 
-0.0007 ±0.0007 


Opw 


n c;r9+0.022 
U.dOZ_o 016 



Table 1. Partial wave contributions to the Olsson sum rule. 

shows that the result is dominated by the contributions 
from S" and P^. The uncertainties in the contributions 
from the partial waves are not symmetric, because the 
input used for the value of the phase shift Sq at 0.8 GeV 
is asymmetric (see appendix IB. 3p . 

Note that the errors are correlated - the uncertainty 
attached to the sum given at the bottom of the table is 
smaller than what is obtained by summing the individual 
errors in quadrature, because it accounts for the correla- 
tions. These also need to be accounted for when evaluating 
the combination 

Opw = Opw - 2a[] ± 5al : (45) 

since the values of the scattering lengths enter the Roy 
solutions, the term Opw depends on these. The net result 
reads _ 

Opw = -0.nOt°oml- (46) 

10.3 Total cross section with It — 1 

In the above notation, the Olsson sum rule amounts to 
the relation Opw ± Or ~ 0. Since the representation ([^5]) 
for Im T'^^'(s, t) is proportional to the residue of the Regge 
pole with the quantum numbers of the p, the cross section 
o''^\s) is proportional to /3p(0) and in addition involves 
the coefficient ri of the pre-asymptotic term, as well as 
the intercept of the p-trajectory, ap(0). Actually, in the 
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narrow range (jl3p we adopt for this intercept, the integral 
is not sensitive to ckp(O). Evaluating it with the central 
value, the result P5|) implies 

/3p(0)(l + 0.104 n) = 80ii^. (47) 
The Olsson sum rule thus leads to a linear correlation be- 
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Fig. 5. Regge term with the quantum numbers of the p. 
The constraints imposed by the Olsson sum rule and by dual- 
ity are compared with the results for /3p(s) obtained from 
other sources: Regge parametrization quoted in the Review 
of Particle Properties [B3] plus factorization; fit to the data 
of Abramowicz et al. [13]; various values taken from the lit- 
erature [2], [37], [Zl]. The ellipse represents the result of our 
analysis in equation l|49)l . 

tween /3p(0) and ri/?p(0), which is indicated by the yellow 
band in Fig. [5] Because of the linearity, the error analy- 
sis simplifies considerably if the variable r\ is replaced by 
ri/3p(0) and we will do that. 

The duality condition, which requires that, on the in- 
terval 1.5 GeV < ^fs < 1.9 GeV, the average over the par- 
tial wave representation agrees with the average over the 
Regge representation yields a second relation of this form: 

/3p(0)(l + 0.351 ri) =6ll^". (48) 

In Fig. [5] this constraint is shown as a blue band. 

The right hand sides of (|T7)) and (1^51) are obtained by 
averaging over the uncertainties attached to the Roy solu- 
tions and are therefore correlated. In order to account for 
these correlations as well as for the uncertainties in ap(0), 
we use the method outlined in section [^ and calculate the 
central values and the correlation matrix of the two quan- 
tities /3p(0) and ri/3p(0). The ellipse shown in Fig. [S] indi- 
cates the corresponding 1 a contour and represents our 
final result for these two quantities. While the bands do 
not account for the correlation between the r.h.s. of (|47l) 
and (|48|) . the ellipse does - this is the main reason why 
the bands are not exactly tangent to the ellipse (the re- 
maining difference stems from the fact that the ellipse also 



accounts for the uncertainty in ap(0), while the bands do 
not). Expressed in terms of the two original parameters 
/?p(0) and ri, our result reads: 

(3,(0) ^Mtll, ri = -0.8tr7. (49) 

There is a modest correlation between these two variables; 
the correlation coefficient is —0.33. The result for /3p(0) 
will be compared with values found in the literature in 
section 110.51 

10.4 Total cross section with It — 2 

The analysis of the total cross section with isospin 2 in 
the t-channel closely parallels the one in section [lO. II The 
grey band in Fig. 4c indicates the sum of the partial waves 
for the corresponding total cross section. The resonances 
are seen clearly, because a*^^^ (s) is not positive and the res- 
onance contributions alternate in sign. The last resonance 
seen is the /93(1690) which generates a small dip. 

The high energy behaviour is parametrized with the 
expression (|?7)) . According to (|^^ . the total cross sec- 
tion falls approximately in inverse proportion to s. In this 
channel, the requirement that, on the interval between 1.5 
and 1.9 GeV, the average over the partial wave representa- 
tion agrees with the average over the Regge representation 
implies 

^g(O) = -24± 29. (50) 

10.5 Discussion: total cross sections 

The Regge parameters emerging from our analysis att — 
are given in equations (1^ . (j^ for It = 0, in (j^ for 
/( = 1, and in ([50| for It — 2. The final result for the total 
cross sections, obtained with the interpolation described 
in section [7] is shown as a solid line in Figs.[T]-4. 

An important feature of our representation is the pres- 
ence of pre-asymptotic terms. As mentioned above, such 
terms become relevant if the Regge representation is used 
at relatively low energies. In the following sections, we 
will show that the semi-local form of duality invoked in 
the present paper does lead to a coherent picture, pro- 
vided the pre-asymptotic contributions are accounted for. 
We emphasize, however, that the simple parametrization 
used for these terms can only give a very rough descrip- 
tion of the contributions occurring underneath the leading 
Regge poles. 

In the cross section with It — 0, the estimate (142]) 
for the pre-asymptotic term reduces the contribution from 
Pomeron and / at 2 GeV by 20± 15%. Fig. 4a shows that, 
towards lower energies, this term widens the band spanned 
by the Regge representation, leaving room for resonances. 
Towards higher energies, the pre-asymptotic contribution 
shrinks. At 5 GeV (lower end of the energy range used for 
the Regge analysis of the irN and NN data in ||33]), it 
amounts to 3 ± 2%. 

Our estimate (|49p for the pre-asymptotic term in the 
channel with /f — 1 is consistent with zero. The main 
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effect of tills term Is to increase the uncertainties in the 
result obtained for the Regge parametrization at low en- 
ergies. Non-leading contributions of this size generate suf- 
ficient broadening, so that the extrapolation below 2 GeV 
becomes consistent with the representation in terms of 
partial waves, at least in the average sense. Note that in 
the channel with It = 1, the resonances generate more 
pronounced fluctuations than for /( = 0, because the con- 
stant background related to the Pomeron is missing here: 
between 1.3 and 1.5 GeV, the sum of the partial wave con- 
tributions drops from 34 ± 2 to 7 ± 3 mb and then again 
increases to 14 ± 4 mb when the energy reaches 1.7 GeV. 

The entire component of the Regge representation with 
It — 2 amounts to a pre-asymptotic contribution. The 
result given in ((50)) shows that, at t = 0, this term vanishes 
within errors. The corresponding total cross section is in 
the range (T^^''(2GeV) — —2 ± 3 mb, small compared to 
the leading cross section cr^^) (2 GeV) = 48 ± 5 mb, which 
is dominated by the contribution from the Pomeron. 

We conclude that the Regge representation exhibits 
the expected qualitative features already at 2 GeV. Our 
results are in general agreement with the total tt^tt^ cross 
section data shown in Fig.[TJ while among the 7r~7r~ data 
above 1.5 GeV shown in Fig.[2j only those of Abramow- 
icz et al. are consistent with our framework. In fact, these 
data provide an independent estimate for the size of the 
coefficient YpTir in the Regge formula (|16p for the to- 
tal TTTT cross section. Taking the values of the coefficients 
B , So, Zj^.^,Yf.^-^ from factorization, neglecting contribu- 
tions with It = 2 as well as pre-asymptotic terms, and 
treating Yp^^ as a free parameter, we obtain an acceptable 
fit to the data of Abramowicz et al. with = 7.7, 
for 5 degrees of freedom and YpTrjr = 19.5 ± 1.6 mb. This 
corresponds to /3p(0) = 100 ± 8, consistent with (H^ . as 
seen also from Fig. [5] 

It can be further checked that the pre-asymptotic terms 
do not distort the description at the energies where the 
data of Abramowicz et al. were taken. Evaluating with 
the results obtained for pi, /?e(0), Pp{0) and ri in 
((50|) , ([49]) , and accounting for the correlations among these 
variables, we obtain — 6.8 for 6 data points, no free 
parameters. As compared to the fit mentioned above, the 
quality thus even improves a little. We conclude that the 
results obtained above for the non-leading Regge contri- 
butions at t = are perfectly consistent with the data of 
Abramowicz et al. 

Figs. 1-4 also show several other Regge parametriza- 
tions of the tttt amplitude at high energies. The red bands 
indicated as PDG 2011 in Figs. 3a-c and 4a, 4b are based 
on the results given in equations (|23|) and (|26)) . which are 
obtained from the irN, NN and NN data above 5 GeV 
with factorization. As discussed in section 14.21 the large 
width of these bands in Figs. 3 and 4b is due to the poor 
extraction of /3p{0). The central value of the factorization 
estimate for /3p(0) is on the high side, but in view of its 
large error bar it is consistent with the value (j49p . ex- 
tracted from totally independent information. 

The Roy equation analysis performed in is based 
on Pennington's representation [5T1. The bands denoted 



as ACGL 2001 in Figs. 3 and 4 show that the Pomeron 
estimate used in [2J, /3p(0) = 46 ± 38, although equipped 
with a large error, is too low. On the other hand, our 
analysis does confirm the value /3p(0) = 72 given in [3T] 
and adopted in 0. As noted in [3], this result is perfectly 
consistent with the theoretical prediction for the combi- 
nation 2aQ — 5aQ of scattering lengths that is relevant for 
the Olsson sum rule. 

The hatched bands that extend down to 1.42 GeV indi- 
cate the Regge representation labeled CFD in [37] , which 
corresponds to (3p{0) = 98.7 ± 1.6, /3/(0) = 31.6 ± 1.9, 
/3p{0) = 58.4 ± 5.5 and l3e{0) = 3.15 ± 7.9. The corre- 
sponding trajectories, as well as the results for the residues 
obtained in the previous works of the same authors, are 
quoted in appendix [X] We note that, although the val- 
ues of the / and p residues are smaller than our results 
in ([^5]) and (HHl) . the Regge representation for the total 
cross section with It — and It = 1 given in (32 is con- 
sistent with ours. The reason is the presence of the pre- 
asymptotic terms in our amplitudes, which are negative 
and in effect reduce the contributions from the dominant 
residues. The larger intercepts of the / and p trajectories 
adopted in [37] also count: at 2 GeV, for instance, with 
the value ap(0) = 0.53 adopted in [37]: the factor s"" is 
larger by 10% than the same factor calculated with our 
value (IT3l) . As concerns the overall behaviour of the cross 
sections above 1.7 GeV, the representation given in [37] is 
consistent with ours. 

Factorization is exploited also in a very recent analy- 
sif[3 of the hadronic cross sections by Halzen et al. |74) . 
which leads to a value Zt^t^ = (12.7 ± 1.4) mb consistent 
with (|17p . but a higher 1/^^ = (16.0 ± 3.9) mb and a 
lower Ypj^T^ — 1.9lJ'Qmb, compared with the numbers in 
(fTS)) . These discrepancies originate in the fact that the au- 
thors of [74] use a smah coupling and a large cou- 
pling Ypjvjv in the denominator of the factorization rela- 
tions (fT5)) . Actually, the values of Yipp and l2pp obtained 
in [74 by independent fits of pp forward scattering data, 
including the total cross section at 7 TeV measured re- 
cently at LHC, are not very different from those given in 
PDG 2011 and quoted in equation But the individual 
contributions of the trajectories are then separated assum- 
ing the equalities = Ya^NN and ypjvjv = Yuinn- The 
values reported in equation ([TU)) and the discussion below 
it show that this assumption is strongly violated. As seen 
from Figs. 4a, 4b and[Sl the unusual Regge picture of the 
TTTT cross sections proposed in [74] is in conflict with our 
analysis. 

Finally, we consider the so-called exchange degener- 
acy, which requires af{t) = ap{t) and (3f{t) = 3/2 f3p{t). 
In the Lovelace-Shapiro- Veneziano model [H] , these rela- 
tions indeed follow from the absence of exotic resonances. 



' Footnote added in August 2012. The comments made 
in this paragraph concern a preliminary version of reference 
[73] , available at http://arxiv.org/abs/1110.1479vl In the 
published version of that paper, which recently appeared in 
Phys. Rev. D85 (2012) 074020, the deficiency is taken care of. 
The new values for the residues V/^-n- = (11.5 ± 0.9)mb and 
VpTTTT = (19 ± 5)mb are consistent with our estimates. 
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In contrast to other analyses, [5] for instance, we do not as- 
sume exchange degeneracy to hold, and it turns out that 
our resuhs for ^/(O) and ,3^(0) given in ^ and dH]), 
respectively, do not have this property. However, as em- 
phasized below equation (^5)) , only the sum of the contri- 
butions from Pomeron and / counts, so that the value 
of /3/(0) is sensitive to the chosen parametrization. To 
illustrate this point, we mention that the parametriza- 
tion of the total NN and ttN cross section adopted by 
PDG before 2002 [58] leads via factorization to the values 
/3p(0) = 57 ± 1, /3/(0) = 111 ± 2, /3p(0) = 82 ± 15, which 
happen to be consistent with exchange degeneracy. We 
conclude that, within our framework, the residue (3f{t) is 
too sensitive to the parametrization used for the Pomeron 
to give an unambiguous answer to the question of whether 
or not the residues of p and / are approximately exchange 
degenerate. 



contributions. For the same reason, however, the specific 
form used for /3f(t) is not crucial. 

The discussion below equation ^ imposes an impor- 
tant constraint on Pf{t). The / is even under charge con- 
jugation, Tf = 1. Hence the absence of ghosts implies 
that /3/(t) must vanish at the place where the trajec- 
tory af{t) passes through zero. This happens around t ~ 
— 0.6GeV^, that is, within the i-interval of interest, close 
to the lower end. The condition leaves little freedom in the 
i-dependence of this residue. We neglect the curvature and 
approximate it with a straight line, /3/(i) oc af{t), so that 
the profile is given by: 



bf = 



a/(0) 



(52) 



Fig. |6] visualizes this approximation: the area spanned 



11 Regge analysis for t 7^ 
11.1 t-dependence of the residues 

The status of the Regge representation for amplitudes at 
nonzero momentum transfer is more uncertain than that 
for the total cross sections. Several global Regge fits of 
data on soft hadronic processes, available in the early lit- 
erature [46], [47], [53] are in rough agreement, but dis- 
agree on details. The absorption corrections, which are 
large for nonforward scattering, are expected to change 
the residues of the pure Regge poles. Since these correc- 
tions are in general process-dependent, they can violate 
factorization [21], [51]. We mention also the more recent 
analysis reported in [71) . However, a systematic treatment 
of all the data sets for i 7^ is not yet available. 

The representations ([21]) . ((25|) . ((27|) leave the i-depen- 
dence of the Regge representation for the imaginary parts 
open. In order to complete our parametrization, we need 
to specify the dependence of the residues Pp{t), /3f{t), 
/3p{t), /3e{t) on the variable t. It is convenient to factor 
out the value at i = 0: 



(51) 



We refer to the functions f3k (t) as profiles. 

In the framework of the Roy equations, we are inter- 
ested in the range -32M^ <t < 4M^, where Mandelstam 
analyticity implies that fixed t dispersion relations hold in 
the interval -AM^ < s < 68M^ ~ (1.15 GeV)^ Numeri- 
cally, the range of interest for the variable t corresponds 
to -0.62 GeV^ < t < 0.08 GeV^. 

In the following, we exploit the fact that the sum rules 
strongly correlate the ^-dependence of the Regge residues 
with the i-dependence of the partial waves. In fact, these 
sum rules can be solved for ^p{t), Pp{t), Pe{t)- The results 
are of the same accuracy as those for the residues at t = 0, 
because the uncertainties in the i-dependence originate in 
the same sources. 

The fourth profile, /3/(t), cannot be determined in this 
way, because the sum rules only concern the sum of Pome- 
ron and / and do not allow us to clearly separate the two 




t (GeV ) 

Fig. 6. Profile of the / residue. 



by the full band corresponds to the uncertainty range in 
([T^ and (|15l) . We can allow for a curvature term - the 
sum rules then yield a somewhat different profile for the 
Pomeron. We have checked that, as long as the curvature 
is not larger than in the other residues, it barely affects our 
results. For comparison, the figure also shows the profile of 
the analogous term in the Regge parametrization of |37j . 
where the symbol P' is used instead of f. 

Concerning the approximation used for the remain- 
ing three profiles, the situation is the same as with the 
parametrization of the phase shifts used when solving the 
Roy equations: since the sum rules allow us to calculate 
these, the parametrization of the profiles is a mere matter 
of convenience. We find that, on the domain of interest to 
us, —0.62 GeV^ < t < 0, a parabolic parametrization is 
adequate: 



/3k{t) = i + bkt + ckr 



P..f,P,e. (53) 



The accuracy to which the sum rules and duality con- 
ditions are obeyed does improve a little if we allow for 
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cubic or quartic terms, but the results can barely be dis- 
tinguished from those obtained with ([53|. In this nota- 
tion, the representation ((5^ for Pf{t) amounts to c/ = 
and the numerical values for intercept and slope of the 
/-trajectory in (fT3)) and dTSl) imply 

6/ = 1.7±0.3 GeV"^ (54) 

As we will see, the central solution for f3p{t) also con- 
tains a zero, not far from the one in 0f(t). This phe- 
nomenon, however, is not related to the absence of a ghost. 
Since the p is odd under charge conjugation, Tp — —1, the 
absence of ghosts calls for a zero in ^p{t) at the value 
of t where ap{t) passes through —1. As this happens near 
t ~ —1.6 GeV^, far outside the range if interest, the condi- 
tion that /3p{t) must vanish there does not impose a signifi- 
cant constraint on our analysis (incidentally, the parabolic 
parametrization ([53| of our central solution does contain 
a second zero at large negative values of t, but it occurs 
outside the region where that parametrization is relevant). 

11.2 Profile of the Pomeron 

The green band in Fig. [7] shows the outcome of our calcu- 
lation for the t-dependence of the Pomeron residue. Nu- 
merically, the upper and lower edges of the band can be 
represented as /3p{t) ± S/3p{t), witf|f| 

^p(i) = 1 + 2.24 1+ 1.76 ^^ (55) 
dj3p(t) = 0.36 1+ 0.26 

For the slope at t = 0, the calculation yields^ 

&P = 2.5±0.4GeV"^ (56) 



|3p(t) 
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t (GeV^) 

Fig. 7. Profile of the Pomeron residue. 



much less rapidly with \t\ than the square of the e.m. form 
factor. 

Results for the structure of the diffraction peak in elas- 
tic nn scattering would be of considerable interest, in par- 
ticular in view of possible direct pion-pion collisions at 
high energies in the near future (for a recent discussion of 
these prospects see [3]). The TOTEM experiment at the 
LHC [92]) which reaches down to very small momentum 
transfer in elastic NN scattering, renewed theoretical in- 
terest in topics like the ratio of the real to the imaginary 
part of the forward scattering amplitude, the comparison 
of the total elastic and inelastic cross sections, the picture 
that represents hadrons as black disks, etc. [T^IMIIM] . 
However, for the reasons given in section HTTl these ques- 
tions are beyond the reach of our framework. 



The figure shows that the sum rule analysis confirms 
the description of the Pomeron profile in terms of the 
e.m. form factor of the pion, 

Mt)=FAtf, (57) 

which was introduced by Donnachie and Landshoff [51] 
and adopted in [60l[6n[67| : the dash-dotted line shows the 
profile obtained from this formula with vector meson dom- 
inance, which predicts F^{t) = 1/(1 - t/M^), M = Mp. 
A parametrization of this form is used also by Cudell et 
al. [7T], but the resulting picture is qualitatively differ- 
ent from the one in f54J: the estimate for the parameter 
M resulting from the fits constructed in [7T] is more than 
twice as large as Mp, so that the Pomeron profile falls off 

* Throughout, the coefficients occurring in the polynomial 
representations of the residues as well as the slope parameters 
bp, bf, bp, be and /3e(0) are given in GeV units. 

^ The value given for bp agrees with P'p{0) ± Sl3p{0) only 
roughly: the former involves the square root of a sum of squares 
of the derivative while the latter is obtained by taking the 
derivative of a square root of a sum of squares. 



11.3 Profile of the p residue 

For the profile of the p, the outcome of our analysis is well 
described by the numerical representation f3p{t) ± 5j3p{t), 
with 

^p(i) = 1 + 3.00 i+ 1.37 (58) 
5Pp{t) = -0.44t + 0.21i2. 

In this case, the value of the slope is given by 

6p = 3.4±0.5 GeV"^ (59) 

Our result for this profile is indicated by the green band in 
Fig.m The residue has a zero at t = -0.42^°;^^ GeV^ The 
uncertainty band partly overlaps with the one describing 
the / residue, which is shown in Fig. [51 Despite the dif- 
ferent origin of the two zeros emphasized at the end of 
section Hi .11 the two profiles are consistent with exchange 
degeneracy. 

The recent Regge analysis of pp and up data in [7T] 
leads to a rather different picture for the f-dependence of 
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Fig. 8. Profile of the p residue, 



the residues. The authors assume factorization, but ne- 
glect the difference between the spin-flip and non-flip am- 
plitudes and do not consider data on the charge exchange 
process. Their representation for /3p(t) contains a zero at 
t = -0.153 ± 0.003 GeV^ or at i = -0.148 ± 0.003 GeV^ 
depending on whether or not a hard Pomeron component 
is included. We recall that a zero at small momentum 
transfer was assumed also in other works, in order to ac- 
count for the so-called cross-over phenomenon, i.e. the 
equality of the 7r+p and ir^p differential cross sections 
around t ~ —0.2 GeV^. The Regge representation of [32] . 
for instance, contains a zero at < = —0.22 ± 0.02 GeV^. 
The analysis of the p-dominated charge-exchange process 
Tr~p — > TT^n in |66) . however, does not confirm the pres- 
ence of such a zero. It appears that the cross-over can 
naturally be explained in the framework of a Regge model 
with strong absorption for ttN scattering [3T1[ST] , without 
the prnr residue passing through zero there. At any rate, 
our results are not consistent with the presence of a zero 
in I3p{t) in the vicinity of the cross-over point. 

Fig. [5] contains three other results for the p profile. 
The one obtained with the LSV model [13] is shown as 
a dash-dotted line. The zero predicted by this model is 
compatible with our result. A similar zero position, at 
— 0.44± 0.05 GeV^, was found in the analysis of Penning- 
ton [21] . which was also based on the Roy equations and 
a set of sum rules (the resulting representation for (3p{t) 
was used as an input in [5]). The figure shows that the 
corresponding band is in very good agreement with the 
outcome of the present analysis. The figure also shows the 
most recent version of the p-profile used by the Madrid- 
Krakow collaboration [37]. Despite some qualitative dif- 
ferences (it is steeper close to the origin and has a larger 
curvature), their result is in fair agreement with ours. 

The uncertainties in /3p{t) rapidly grow with \t\: note 
however that what counts in the applications is not Pp{t), 
but the imaginary part of the amplitude which is propor- 
tional to P p{t){s / si) "p(*). At an energy of 2 GeV, the term 
(s/si)"''^*) drops roughly by a factor of 2 as t runs from 



the upper to the lower end of the interval relevant for the 
partial wave projections, > t>— 0.6 GeV^. 



11.4 Effective residue in the channel with It — 2 

Within errors, the effective residue in the exotic channel, 
(3e{t), vanishes at i = 0. The corresponding profile Pe{t) 
does therefore not represent a meaningful quantity. In- 
stead, Fig.jS] shows the t-dependence of the residue itself. 
The upper and lower edges of the band shown can be rep- 
resented as (3e{t) ± 6/3e{t), with 

/3e(i) = "24- 139^-94^2, (60) 
SPeit) = 29 -2t- 19^2, 

and for the slope at t = 0, we find 

/3^(0) = -120±40 GeV"^ (61) 

In the entire range considered here, this residue is small. 
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Fig. 9. t-dependence of the effective residue with It = 2 



The uncertainty is approximately independent of t. The 
value at i = is related to the size of the total cross section 
in the channel with It — 2, which is shown in Fig. 4c. At 2 
GeV, this cross section amounts to a^^\s) = —2 ± 3 mb, 
small compared to the leading term from the Pomeron. 
At negative values of t, the residue passes through zero 
and then becomes positive. In view of the fact that the 
amplitude is proportional to /3e(i) s"^ the uncertainty 
even decreases a little towards the lower end of the range 
of t- values where our analysis applies. 

Fig.ini also shows the residue used in the analysis of 
[3S1I37], where the t-dependence of the residues is taken 
from [221117]. In the latter references, the Regge parame- 
ters for TTTT scattering are determined on the basis of fac- 
torization. There is, however, barely any direct informa- 
tion on the component of the Regge representation with 
It = 2. For what concerns /3e{t), the expectation expressed 
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in [3S1I37] that their Regge parametrization describes ex- 
perimental data in the region 1.42 GeV < y/s < 20 GeV, 
—0.4 GeV < t < AM^ merely represents an educated 
guess. Our quantitative analysis does confirm that I3e{t) is 
small, but it does not support the uncertainties attached 
to their parametrization, which in our opinion are grossly 
underestimated. 



In Fig. llOl the low energy contributions to the function 
S{t) defined in (p9| are shown as a grey band (numerical 
values in GeV units). To check whether the high energy 
contributions do approximately cancel these, we flip their 
sign, so that instead of canceling, the two contributions 
would become equal if the sum rule was obeyed exactly. 
The result is indicated by a hatched band. Evidently, the 
sum rule is very well obeyed within errors. 



12 Consistency checks 

The present section concerns the internal consistency of 
our results: we compare the contributions to the sum rule 
and duality relations arising from the low energy region, 
where we are using the partial wave representation, with 
those from high energies, for which we use the Regge repre- 
sentation. For the relations to be approximately satisfied, 
the two contributions must approximately cancel. 



12.1 Checking the Olsson sum rule 

Consider first the Olsson sum rule. The contribution from 
the integral over the partial waves and from the subtrac- 
tion constants was analyzed in section [I0.21 with the result 



O 



PW 



-0.110 



+0.022 
-0.015- 



Evaluating the remainder with the 



Regge represention constructed in the preceding sections. 



we obtain On = 0.106 



-fO.016 
-0.018- 



The two terms thus indeed 



cancel within errors. This does not come as a surprise, be- 
cause the Olsson sum rule occurs among the conditions im- 
posed when constructing the Regge representation. What 
the agreement mainly checks is that the Olsson sum rule 
did receive enough weight among the many constraints 
considered, so that the outcome for the Regge representa- 
tion does respect it. 



12.2 Verifying the sum rule S{t) — 




-0.3 -0.2 

t (GeV^^ 



Fig. 10. Validity of the sum rule S{t) = 0. 



12.3 Checking crossing symmetry 

As discussed in section |S1 the fact that the crossed chan- 
nels of elastic nn scattering are identical imposes strong 
constraints on the scattering amplitude. In particular, these 
imply a set of sum rules of the form Ck{t) = 0, where 
k = 0,1,2 and t is a free parameter (see appendix [C|) . 
In Fig. [Til the red and blue bands for Co(i) and Ci{t) 
show the partial wave contributions to the integrals oc- 
curring in the definition (j75p of these functions, while the 
hatched ones indicate the corresponding Regge contribu- 
tions, changed in sign (numerical values in GeV units). In 
order to disentangle the bands for Ci{t) and C2{t), the 
sign of the low energy contributions to C2{t) is flipped 
instead of the one from the Regge region. The fact that 
the partial wave and Regge bands are consistent with one 
another demonstrates that the sum rules Ck{t) = are 
well obeyed within errors. 




t (GeV ) 

Fig. 11. Validity of the crossing symmetry sum rules. 



12.4 Checking duality 

In the construction of our Regge representation, we made 
use of a semi-local form of duality: we required the average 
over the interval 1.5 GeV < ^/s < 1.9 GeV of the partial 
wave representation for the total cross sections to agree 
with the average of the Regge representation over the same 
interval. Moreover, we imposed the same condition also at 
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nonzero values of t, replacing the cross sections by the 
quantities p(s) ImT^^*) (s, t). In Fig.[T2l the fuh bands in- 
dicate the averages over the partial wave representations 
for these quantities, while the hatched bands represent the 
averages over the Regge representations (numerical values 
in mb). The figure shows that the two averages agree: 
our partial wave and Regge representations do satisfy the 
semi-local form of duality invoked in the present paper. In 
fact, since we required the duality conditions for the total 
cross sections to be strictly obeyed, the full and hatched 
bands strictly agree at t = 0. 



13 Summary and conclusions 

In the present paper we have analyzed the high-energy tttt 
scattering amplitude. Our main motivation is not high en- 
ergy per se but rather that this part of the amplitude plays 
a role in the Roy equations, even when these are used to 
determine the amplitude at low energy. Close to threshold 
the relevance of high energy is rather limited, but as one 
moves up, towards 1 GeV and higher, a precise determina- 
tion of the high-energy contribution becomes important. 

We have adopted a Regge parametrization to describe 
the TTTT scattering amplitude at high energy. Some of its pa- 
rameters are well known like, e.g. the Regge trajectories. 
Moreover, some Regge residues can be determined from 
high-energy ttN and -A^-A^ scattering data [6TI - l63l[7T] as- 
suming factorization. Several other parameters, however, 
are poorly or not at all determined. Aim of the present 
paper was to pin these down. In order to do this we have 
worked under the assumption that the partial wave and 
Regge representations can be joined smoothly in the re- 
gion between 1.7 and 2 GeV - not only a.t t — but also 
as a function of t, down to about —0.6 GeV^. In this set- 
ting we have imposed a set of sum rules which follow from 
general principles, like crossing symmetry. 

While 2 GeV can hardly be considered asymptotic, it 
has been already observed in the literature (see [61] and 
references therein) that, at least on average, duality be- 
tween the partial wave and Regge representations works 



well even at energies that low. In order to both partly test 
this assumption and make our results less dependent on 
it, we have supplemented the standard Regge parametriza- 
tion with pre- asymptotic terms (terms that are suppressed 
at asymptotic energies with respect to the dominating 
Regge contributions) and have determined them by re- 
quiring that a smooth transition between the partial wave 
and Regge representations does take place. 

For the partial wave representation, we use a simple 
parametrization which approximately describes the data 
between 1.15 and 2 GeV. Below 1.15 GeV, we solve the 
Roy equations - in this region, the exact form of the 
parametrization becomes irrelevant. In order to solve the 
Roy equations, we need an input for the high-energy part, 
which, on the other hand, is what we aim to determine. In 
summary, we have simultaneously solved a number of con- 
straints - Roy equations, duality and sum rules - which 
impose relations among the different partial waves at low 
energy and the high-energy amplitude. This is not an easy 
task, and to perform it in a way which is as practical as 
possible we have adopted an iterative scheme: 

1. take an input for the Regge parameters; 

2. solve the Roy equations for all partial waves up to and 
including the G waves; 

3. impose the duality and crossing symmetry constraints 
in order to determine the Regge parameters; 

4. if the change in the Regge parameters is significant, go 
back to 1. 

The outcome of this analysis is an improved Regge para- 
metrization of the TTTT scattering amplitude at high energy, 
which: 

1. is consistent with the high-energy data on hadron- 
hadron scattering under the factorization hypothesis; 

2. satisfies the Olsson sum rule and its extension to i 7^ 0, 
as well as a set of other sum rules which follow from 
crossing symmetry; 

3. smoothly joins the partial wave representation in the 
region just below 2 GeV, also as far as the t dependence 
is concerned. 

A summary of the results of the present analysis and 
of the parameters taken as input is presented in Tabled] 
Besides the constants /3p(0),pi, ri, /3e(0), which are rele- 
vant for the total cross sections, our analysis also provides 
the t-dependence of the residues, which we parametrize 
in terms of quadratic polynomials. We do not repeat the 
explicit results for the profiles Px{t) = Px{t) I Px{^) here. 
They can be found in the equations listed in the table. We 
stress that our method for determining the t-dependence 
of the residues is new and puts their representation on a 
more solid basis. 

The main novelty of our analysis is twofold: First, 
the introduction of pre-asymptotic terms, which we have 
shown to be needed in order to reach a smooth transi- 
tion between the Regge and partial wave representations, 
especially in the channel with It = 0. Second, the determi- 
nation of the p residue via the Olsson sum rule, which is 
much more robust than factorization. These two elements 
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explain some of the differences found with earlier analy- 
ses. For a detailed comparison of our results for the total 
cross sections, both with data and the literature, we refer 
the reader to section [TITSl 

As we have discussed in section [11.21 the outcome of 
our calculation for the Pomeron profile is consistent with 
some of the results available in the literature, but differs 
from others. Within the accuracy of our calculation, this 
profile is well described by the square of the pion electro- 
magnetic form factor. 

Our results for f3p{t) indicate that this profile has a 

zero in the range t = — 0.42^q'24 This confirms the 

results of Pennington JTD , which in fact were based on a 
similar approach. On the other hand, as discussed in sec- 
tion 111.31 the outcome of our calculation is not consistent 
with parametrizations that explain the cross-over in ttN 
scattering in terms of a zero in /3p{t). 

As shown in section 111.41 the sum rules studied in 
the present paper also provide a good handle on the t- 
dependence of the Regge amplitude with It — 2. The main 
result here is that, on the interval > t>— 0.6 GeV^, this 
amplitude is small compared to the contribution from the 
Pomeron. The size is comparable to the pre-asymptotic 
terms in the non-exotic channels with /( = 0, 1. As is the 
case with these, the absolute size of the amplitude with 
/( = 2 is subject to large uncertainties. 

New solutions of the Roy equations for the S, P, D, F 
and G waves based on this high-energy input have been de- 



Parameter 


Value 


Source 


ap(0) 


1 


PDG 2011 


«/(0) 


0.54 ±0.05 


our estimate 


«p(0) 


0.45 ±0.02 


our estimate 


a40) 





our estimate 


a'p(0) 


0.25 ±0.05 GeV"^ 


diffraction peak 


a'fiO) 


0.90 ± 0.05 GeV"^ 


Chew-Frautschi plot 


«p(0) 


0.91 ±0.02 GeV"^ 


Chew-Frautschi plot 


c^'M 


0.5 ±0.1 GeV-2 


our estimate 


B 


0.025 ± 0.001 


PDG 2011 


\/so 


5.38 ± 0.50 GcV 


PDG 2011 


Pp{0) 


94 ± 1 


factorization 


/3/(0) 


69 ±2 


factorization 


Pfit) 


equation (1521) 


linearity, no ghosts 


/3p(0) 




Olsson s.r., duality 




-24 ± 29 


duality 


Pi 


-0.8 ±0.6 


duality 


ri 


-o.8t°:? 


Olsson s.r., duality 


bp 


2.5 ±0.4 GeV-2 


sum rules, duality 


bp 


3.4 ±0.5 GeV"^ 


sum rules, duality 


P'M 


-120 ±40 GeV-2 


sum rules, duality 


Mt) 


equation (1551) 


sum rules, duality 


Mt) 


equation (1581) 


sum rules, duality 


P4t) 


equation (1601) 


sum rules, duality 



Table 2. Regge representation of the scattering amplitude. 
The upper part lists the range used for the input parameters, 
while the lower part indicates the results of our calculation. 



rived and will be presented in a forthcoming publication. 
These solutions extend up to 1.15 GeV and are suitable to 
be used in dispersive analyses of other quantities and pro- 
cesses, like the hadronic vacuum polarization contribution 
to {g — 2)^, rj — Stt, Kaon decays in multipion final states 
etc. 

These extended Roy solutions also provide an improved 
basis for the determination of the pole positions of reso- 
nances decaying into two pions, like the a or the p. We em- 
phasize that extrapolations based on explicit parametriza- 
tions of the energy dependence, such as those in [^S] or 
in |96j . cannot compete with the analytic continuation 
provided by dispersion theory. In the case of resonances 
with a large width, the result for the pole position is sub- 
ject to large theoretical uncertainties associated with the 
choice of the parametrization, even if those on the real axis 
are small [ST]- The result for the mass and width of the 
(7 obtained by extrapolating the explicit parametrization 
in [95], for instance, differs significantly from the pole po- 
sition calculated on the basis of the Roy equations [3] , de- 
spite the excellent agreement on the real axis. The disper- 
sive method has the advantage that it provides a complete 
error analysis which includes, in particular, also the uncer- 
tainties associated with the presence of inelastic channels. 
The calculations performed in [Ml [22] corroborate this 
statement: the result for the pole position agrees with [3], 
within the quoted errors. 
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A Literature on the Regge parameters 

Most of the literature on the Regge parameters deals with 
NN and ttN scattering, for which fits of various experi- 
mental data sets were performed. The parameters relevant 
for TTTT scattering are in general obtained only indirectly, 
using factorization [47,65i7ni74|. Several modifications of 
the standard expression ((!]) were also considered in the lit- 
erature: for instance, an additional factor r{a(t)) is often 
introduced in the denominator of ([Ij in order to remove 
the unphysical poles due to sinTra. Other factors depen- 
dent on the trajectory, which remove unphysical poles, are 
adopted in [17], [21I3113I], [5a[7S], etc. 

In the present review, we mainly consider the trajec- 
tory parameters, the residues at i = and the slope of 
the Pomeron profile at this point. For simplicity, we mostly 
quote only the central values, omitting the errors. We men- 
tion that, although the values reported in the literature 
for a certain parameter may be quite different, correla- 
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tions with different values of other parameters may lead 
to a similar overall description. 

A.l Pomeron 

In the early literature, the contribution of the Pomeron 
to the total cross sections was assumed to be energy inde- 
pendent, which in view of ([5]) corresponds to a trajectory 
with intercept ap(0) = 1. In equation ([5]), the Pomeron is 
described by a triple Regge pole with ap(0) = 1, leading 
to a cross section that contains a constant term and a term 
which grows with In^s. The "soft Pomeron" [5 5) with an in- 
tercept slightly larger than 1 (for instance, ap{0) — 1.093 
[58] or ap(0) = 1.08 [61]), also yields a very good descrip- 
tion over a very wide range of energies, while a "hard" 
Pomeron with ap(0) = 1.45 was shown [57j to explain the 
data on ep scattering at HERA. In [71] both of these vari- 
ants of the Pomeron are used for the description of pp, np 
and Kp elastic scattering at small t. 

Information on the t-dependence of the trajectory is 
obtained from the slope of the diffraction peak in pp and 
TTp collisions at high energies: 

Bap{s)^dt\n[da-P{s,t)/dt]^^„, (62) 

where a stands for p,p or tt^. The quantity (j62p is ex- 
pressed as 

Bap{s) = 2[b''/ + a'p(O) ln(s/si)] , (63) 

in terms of the slope of the trajectory, a'p{0), and the 
slope 6p'' of the Pomeron profile at t = 0, defined in the 
case of TTTT scattering in equation (j53p . 

The slope of the diffraction peak in proton-proton scat- 
tering was measured by many experiments at various en- 
ergies (for references on the data sets see [7T1[72]). Re- 
cently, the slope Bpp was measured at ^/s = 7 TeV by the 
TOTEM Collaboration at LHC [S2]- The trajectory slope 
ap(0) = 0.25 was shown to be consistent with high energy 
data, in particular in the context of the soft Pomeron [SI] . 
A recent comprehensive analysis 1721 , based on the avail- 
able experimental data on pp and pp scattering over a large 
energy scale, reports the average value ap(0) = 0.24. 

Recent choices of the Pomeron trajectory are: 

ap{t)= 1.0-h0.20i [37], 1.09 + 0.33 1 [71], 1.08 -h 0.25 1 [75], 

to be compared with the input used in the present paper, 
ap{t) = 1 -h (0.25 ± 0.05) t. 

The Pomeron residue /3p(t) = (3p^{t) is known with 
less precision. For the value at t = 0, we quote: /3p(0) = 
46 [g, 118 [32], 99 [37], 66 [65], 98 [71], to be com- 
pared with the input we are using, /3p(0) = 94 ± 1. Fig. 
4a. illustrates the fact that the value used in ^ is on 
the low side, while the one in |74j is high. The value de- 
rived in [BS] (which is based on factorization and the older 
parametrization of NN and ttN data given in [SS]) is also 
low, but this is partly compensated by a higher value for 
/3/(0). 

As concerns the t-dependence near t — 0, an exponen- 
tial form 

(3f{t)^l3f{0)e'''^^'\ (64) 



was often assumed for various types of elastic reactions 
ab — ab (early references are f45','47','50^ , for more recent 
work see for instance 5^32,34-37 , [66,75 ). For the slope in 
TTTT scattering, bp = &p'^, factorization implies 

bp - 267 - ^p ■ (65) 

Since the 7rp data at t ^ are less accurate, the extrac- 
tion of 6p from this relation is not very precise. The val- 
ues proposed in the literature scatter wildly: 6p = 4 [TO] . 
3.3 [37], 0.42 [13, 0.79 [71], to be compared with the 
result of our calculation, bp = 2.5 ± 0.4. 



A. 2 / pole 

The parametrization of the / Regge pole (denoted also as 
P' in the literature [3^[31ti371H7] ) is linked with the one 
of the Pomeron. Indeed, both contribute to the amplitude 
with It = and are largely indiscernible. The relative 
weight of the contributions from Pomeron and / depends 
on the parametrization chosen for the Pomeron. For one 
of the fits obtained by Rarita et al. 23, for instance, the 
Pomeron dominates over the / already below 1.8 GeV, 
while for the others this is the case only at much higher 
energies. 

Recent parametrizations of the /-trajectory are: 

af{t) = 0.53-h0.90t [37], 0.69-h0.80t [SS], 0.96 -h 0.58 t-h 
0.03 1^ [5i, 0.61 + 0.82t [71], 0.71 + 0.83t [75], while we 
are using af{t) = (0.54±0.05)-h (0.90±0.05) t as an input. 

The values of /3/(0) proposed in the earlier literature 
[^H7] cover a large range. Some values adopted recently 
are (compare section dnS]) : /?/(0) 108 [5], 32 [37], 103 
[55] . 153 [7T], 123 [71], to be compared with the input we 
are relying on, /3/(0) = 69 ± 2. 

A. 3 p pole 

The p trajectory is relatively well known. Some choices 
made in the literature are: 

a p{t) ^ 0.53+0. 90 1-0. 15 1'^ [37], 0.44-^0.96^-0.02*^ [Sg, 
0.48-H0.88t [59], 0.53-F0.8f [Ml|75], 0.47-h0.91f [71], 
while our estimate reads ap{t) = (0.45 ± 0.02) + (0.91 ± 
0.02) t. 

For the residue, the values /3p(0) « 104.2 given by 
Kane and Seidl [FT], /3p(0) w 107.6 of Irving and Worden 
[52] and /3p(0) w 116 predicted by the Lovelace-Shapiro- 
Veneziano model [2] are on the high side. Pennington [3T] 
obtains /3p(0) « 72, a value which was adopted also in [5]. 
These numbers are to be compared with the outcome of 
our calculation, which reads /3p(0) — 84li[g. 

A smaller value, /3p(0) = 33.2 was proposed in [35], 
invoking consistency with data on tttt scattering at high 
energies. In the subsequent papers [31U37| . the authors 
resorted to dispersive analyses of the pion-pion amplitude 
and gradually increased their estimate: /3p(0) = 37.1 [34| . 
48.2 [35], 58.4 [37]-CFD solution. The difference between 
the values obtained in [31]- [3B] and our result is due partly 
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to a different shape of the isoscalar S wave below 1 GeV 
(the phase shift Sq{s) given in these works displays an 
unphysical "hump" around 0.8 GeV, see the discussion 
in [5711100] V We mention that in the most recent result of 
this group, quoted as CFD solution in [37], this hump is 
no longer present. 

Other recent values are /3p(0) — 84 (obtained in [65] 
with factorization, using the older parametrization of cross 
sections in [5S]), which is consistent with our result, and 
the very small value /3p(0) = 9.8 found in [74], which is 
discussed in section [TU31 

The results given in the literature for the residue /3p (t) 
at t are rather controversial. The values obtained from 
dispersion sum rules for the tttt amplitude in the seven- 
ties [T71]18II21| were not conclusive, due to the poor quality 
of the data at low energies. The phenomenological infor- 
mation provided by factorization was also rather poor: 
since the contribution of the p-pole to the forward NN 
amplitudes is small, it was neglected altogether in some of 
the early fits [Uj . The studies of the p-dominated charge 
exchange process tt~p tt^ n indicated that indeed the 
spin-flip pNN coupling is larger than the non-flip cou- 
pling, which is the only one contributing to the forward 
NN scattering amplitude. This feature is confirmed by 
the recent analysis of the pion-nucleon charge exchange 
amplitudes [55] . 



the sum rules require a zero and a comparatively large, 
negative slope, P'^{0) = -120(40). 



B Roy equations 

B.l Fixed-t dispersion relation 

The derivation of the Roy equations (see [T] and appendix 
A of [2]) is based on the fixed-f dispersion relation obeyed 
by the s-channel isospin amplitudes T = {T'^ ^T"^) in 
the interval -28M^ <t< 4M^, written as [T] 

Tis,t) =. (4M2)-i(sl + iC,t+uC,„)r(4M2,0) 

+ / ds' g2{s,t,s')lmT{s',0) (66) 
Jam?. 



/ ds'g3(s,i,s')IinT(s',t) . 

J AMI 



The subtraction term is fixed by the S'-wave scattering 
lengths: 



T(4Af2,0)=32^(aS],0,a2), 



the crossing matrices Cm 
are given by 



Cut J C — Cus^ Cs 



(67) 



Cu 



A. 4 Regge cuts 

The Regge cuts are introduced in phenomenological anal- 
yses especially for extending the validity of the Regge 
model to larger values oi \t\ (see for instance the recent 
works [551IBB] ) . As we are interested in a rather small range 
of we consider a Regge cut only for the description of 
the exotic amplitude with It — 2, where it is the dominant 
contribution. 

The information about this term is poor, since it does 
not contribute above 5 GeV, where the Regge fits to the 
NN and ttN amplitudes are usually done ^63. 71. 73i.74j. 
Alternative treatments based on sum rules for tttt scatter- 
ing in the seventies were not conclusive. Pennington [21] 
notices a great sensitivity of the finite-energy sum rule 
for the amplitude with It = 2 with respect to the transi- 
tion point from low to high energy regimes. He assumes 
/3e(0) « 0, this value being adopted also in [2]. A small 
effective residue Pe{t) appearing in equation (1271) can oc- 
cur from cancellations of the p — p and 02 — 02 cuts, as 
argued by Worden [48] , using symmetry arguments in the 
Reggeon-calculus of Gribov [49]. In [37] the authors ob- 
tain ^e(O) = 8 in the UFD fit, and /3e(0) = 3 in the CFD 
fit. This is consistent with our result, which comes with a 
large error: /3e(0) = -24± 29. 

The t-dependence of /3e(0 is also poorly known. In 
[37], the residue is taken to be proportional exp&et, with 
he = 2.4. The ansatz prevents Pe{t) from passing through 
zero and makes it shrink if the momentum transfer grows 
(compare Fig.|9]). The numerical value for the slope is 
/3^(0) = 20 and ^^(0) = 8 for the UFD and CFD fits, 
respectively. Our framework does not support this ansatz: 



Cu 



C'st 



C su 



and the kernels g2{s,t,s'), 53(5,^,5') represent 3x3 ma- 
trices built with Cst, Ctu and Cs„, 



92{s,t,s') = - 



93{s,t,s') = 



t 



TTS' (s'-4Af2) 
1 , Csu 



{uCst + sCstCtu) X 



S' — t S' ~ Uq ^ 

su / 1 



Csu 



7rs'(s' — Uq) \s' — s s' — u J ' 



.f2 



where u — 4M^ ~ s — t and uq — 4M, 

We add a comment concerning the subtractions. As 
emphasized in [37], a single subtraction suffices: the 01s- 
son sum rule (137)) can be used to express one of the two 
subtraction constants in terms of the other. The resulting 
set of relations is referred to as the GKPY-equations. Con- 
versely, the dispersive representation for T{s,t) obtained 
in that framework [37] differs from the one in (|66p by a 
multiple of the Olsson sum rule. The advantage of using 
a single subtraction is that this allows a better determi- 
nation of Oq and Aq from experiment (see the discussion 
at the end of section [5]). The price to pay is that the dis- 
persive representation for the partial waves then becomes 
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more sensitive to tlie high energy part of the input used, 
because the GKPY-equations converge less rapidly than 
the Roy equations. As we are relying on the theoretical 
predictions for Qq and Oq, retreating to a single subtrac- 
tion is no gain for us. 

B.2 Dispersive representation of the partial waves 

Exploiting crossing symmetry, the partial wave projection 
of the amplitude in (|5T|) can be written in a form that 
involves a smaller range of t-values: 

tl{s) = ^ f dzP,{z)T\sM)- (69) 

This step extends the range of validity of Roy's represen- 
tation [1], 

2 oo /.oo 

t\{s) = klis) + ds' K'.fis, s') Im4'(s') ■ 

(70) 

The term kj{s) denotes the partial wave projection of the 
subtraction term, which shows up only in the S- and P- 
waves, 

fci(s) = ai 5° + '-^^ (2a° - 5al) x (71) 

The kernels X// (s, s') are explicitly known functions (see 
appendix A of jl]). They contain a diagonal, singular Cau- 
chy kernel that generates the right hand cut in the par- 
tial wave amplitudes, as well as a logarithmically singular 
piece that accounts for the left hand cut. The validity of 
these equations has rigorously been established on the in- 
terval - 4M2 <s< 68Af2 = (1.15 GeV)2. 

In [2], Roy's equations were solved for the S- and P- 
waves in the region ^/s < 0.8 GeV. In the present paper we 
solve these equations in the whole range -y/s < 1.15 GeV, 
including also the higher partial waves, which become in- 
creasingly important as the energy grows, because they 
are given more weight in the sum and the angular mo- 
mentum barrier becomes less effective: for instance, above 
1.6 GeV, the S'-wave accounts for less than a quarter of 
the total cross section with = 0. In order to have a 
reliable description, we solved the Roy equations for the 
waves with £ < 4 below 1.15 GeV and made use of the fact 
that, below 2 GeV, the masses, widths and tttt branching 
fractions of the resonances are now known quite well 6 I . 
A detailed description of the solution of the Roy equa- 
tions will be given elsewhere [35] ~ in the following, we 
limit ourselves to a cursory discussion. 

B.3 Input used when solving the Roy equations 

The input required for the construction of the solutions of 
the Roy equations was outlined in section[5] There, we also 



reviewed the status of our knowledge of the subtraction 
constants. In the present appendix, we briefly discuss the 
other elements of the input used for the partial waves. 



B.3.1 Phases 

Our analysis relies on the following phenomenological es- 
timates for values of the S- and P-wave phase shifts in the 
low energy region (^ij = 0.8 GeV, y^s^ = 1.15 GeV): 

6°{sa) = 82.3°il°° , (72) 
S^iAM]^) = 185° ± 10° , 
<5!i(w) = 260°±10°, 

SI{sa) = 108.9° ±2°, 
SliSma.) = 166.5° ±2°. 

We briefly comment on these estimates. 

The phase shift analyses available at the time were 
used in [5] to estimate the S° phase shift at 0.8 GeV, with 
the result 6q{sa) = 82.3° ± 3.4°. This estimate was criti- 
cized in |32| . however, and significantly higher values were 
suggested. We did reply to this objection [Sj, but when 
analyzing the consequences of the representation for the 
S°-wave for the mass and width of the a [39^, we decided 
to inflate the uncertainty range by a factor of 2, in order 
to make it evident that the outcome is not tied to an input 
which at the time appeared to be controversial. 

In the meantime, the dust settled - the discrepancy 
disappeared: the result Sq{sa) = 85.7° ± 1.6° obtained 
in the most recent paper of the Madrid-Krakow collab- 
oration (Table XII in |37j) differs from the old estimate 
in 2 by less than one standard deviation. Moreover, the 
thorough analysis of the experimental situation concern- 
ing the phase shift in [95] fully confirms the results 
of [1113], also with regard to the value at 0.8 GeV. Last, not 
least, we point to a recent paper by Moussallam [52], who 
has analyzed the problem, solving the Roy equations with 
matching point at 2Mk and relying on the experimental 
results for the scattering lengths obtained by NA48/2 [83] 
to fix the subtraction constants. He shows that, currently, 
the uncertainties in the value of 6q{sa) are dominated by 
those in the input used for the elasticity above KK thresh- 
old: while the shallow dip in 77§(s) obtained from the data 
on the inelastic channels leads to Sq{sa) = 82.9°±1.7°, the 
deep dip indicated by the data on elastic scattering yields 
S^{sa) = 80.9° ± 1.4°. Both of these resuhs are within the 
range estimated in 2,. Hence we might just as well return 
to that estimate, but, for the time being, we prefer to stick 
to the broad uncertainty range adopted in [39J - we plan 
to examine the present experimental situation concerning 
the low partial waves in a forthcoming article 38 . 

Phase shift analysis indicates that 6q{s) passes through 
180° in the vicinity of KK threshold, but it is still not 
known for certain whether this happens below or above 
that point, although the second option appears to be more 
likely. The error bar attached to our estimate allows for 
both possibilities. The values for 5q(4M^) given in [55] 
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are consistent with our estimate in (j72p . albeit somewhat 
higher. 

The P^-wave is experimentally very well determined 
through the data on the e.m. form factor of the pion, so 
that our estimates for the phase shifts for SI{sa) and 
^i(smax) in ((7^ come with small errors. 

The list ([7^ does not contain entries for the exotic 
phase shift Sq{s). The reason is the following. As discussed 
in [2], the solutions of the Roy equations in general de- 
velop a cusp at the matching point, which in the present 
framework is taken at ^fs^^^^ =1.15 GeV. We remove the 
unphysical phenomenon by requiring the phase shifts as 
well as their derivatives to be continuous at the matching 
point. This in effect imposes a condition on the behaviour 
of the imaginary parts above this point: we find acceptable 
solutions only if we do not prescribe the value of (5g (smax), 
but let it float (for a more thorough discussion, see [5]). 
The same holds for the higher partial waves: the require- 
ment that the phase shift and its derivative are continuous 
at the matching point admits solutions only if the value 
of the phase at this point is not prescribed. 

For the numerical evaluation of the Roy equations, we 
need an explicit parametrization of the phases. Since the 
solution is uniquely determined, the choice of the parame- 
trization is a matter of convenience - the only requirement 
to be met is that it is flexible enough to come close to 
the exact solution. In the case of the S°-wave, the fact 
that there is a resonance in the immediate vicinity of KK 
threshold entails a rather intricate behaviour there and 
special care is needed to find solutions of good quality, but 
for the other waves, we did not encounter such problems. 
For details, we refer to |38) . 



a Breit-Wigner parametrization and calculate the corre- 
sponding inelasticities from mass, width and branching 
fraction for the decay into tttt quoted in the data tables. 
The exotic waves and are driven by the other par- 
tial waves. Since their inelasticities are very small, these 
barely affect our results. For a more thorough discussion, 
we refer to l38l. 



B.3.3 Imaginary parts 

Although phenomenology does reduce the range allowed 
by unitarity, < Imt^(s) < l/p{s), the uncertainties in 
the behaviour of the partial wave imaginary parts above 
the matching point are considerable. 

We use a crude parametrization involving Breit-Wigner 
terms for the resonances and polynomials for the back- 
ground, impose continuity of function and derivative at 
the matching point and allow for one free parameter in 
the representation used, for each one of the partial waves 
with £ < 4. It is convenient to identify the free parame- 
ter with the value of the corresponding imaginary part at 
y/Jb = 2 GeV. Our estimates for the latter read: 

Imi[](s6) = 0.4 ±0.25, Imi}(sfc) = 0.26 ± 0.13 , (73) 

liatlisb) = 0.3 ± 0.15 , liat%Sb) = 0.2 ± 0.075 , 

ImiJ(sh) = 0.23 ± 0.07 , Imt^(sb) = 0.12 ± 0.06 , 

ImtKsb) = 0.22 ± 0.03 , ImtKsb) = 0.01 ± 0.01 . 

We plan to describe the parametrization used in more de- 
tail [38;. 



B.3.2 Elasticities 



C Sum rules from crossing symmetry 



The input used for the elasticity of the partial waves in 
the region where the Roy equations are solved plays a sig- 
nificant role, not only in resolving the structure of the S'^ 
wave in the vicinity of KK threshold, for instance, but 
also when analyzing the behaviour of the partial waves 
off the real axis, in order to determine the position of the 
resonance poles, such as the one of the a. The elasticities 
ensure that the effects on the elastic tttt scattering ampli- 
tude generated by inelastic processes like tttt ^ KK, for 
instance, are properly accounted for. 

The interference between the KK cut and the pole 
associated with the /o(980) generates a dip in the elas- 
ticity of S", immediately above ^/s = 2Mk- In the region 
between 2Mk and 1.15 GeV, where it matters for our ana- 
lysis, the uncertainties in the elasticity of S° are gradually 
becoming smaller. For a recent, very thorough analysis of 
this field, we refer to [22] ■ We describe the behaviour by 
means of the Flatte formula flOl.l and vary the parame- 
ters contained therein in a rather broad range, so that the 
uncertainty band for ?7g(s) generously covers the available 
parametrizations. The inelasticities of and S^ are also 
estimated on the basis of available phase shift analyses. 
For the higher partial waves with Is = or 1, we use 



A set of sum rules for fixed t is obtained by taking the 
first derivative with respect to s of the l.h.s. of (pij) and 
then setting s = 0. Since the resulting expression has a 
zero both at t = and at t = 4M^, it is convenient to 
divide by the factor t{t - 4:M^). The sum rules then take 
the form 

Ck{t)^0, fc = 0,l,2, (74) 

where Cq (t) , Ci {t) , C2 (t) are integrals over the imaginary 
part of the scattering amplitude: 



Ck{t) 



4M2 
00 



dsnk{s,t) lmT^''\s,t) 

nOC 

/ dsfk{s,t)lmT\s,0) 

dshki{s,t) dt'lmT^ {s,t')\ 



(75) 



2 

E 

i=0 



4A/2 



The functions ImT ^°\s,t), ImT(i)(s,t), lmT^'^'>{s,t) are 
defined in (|40)) . The first line shows that the sum rules 
decouple, in the sense that the i-dependent part of the 
scattering amplitude enters with definite isospin in the 
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t-channel. The explicit expressions for the corresponding 
kernels read 



2s 4 


- 1 - AMI 




f t - AMIY 


t 


- AMI 


TT s'^{s - 


f t - AM-^f 




' t - AMI 




f t - AMif 



no{s,t) 
ni{s,t) 
n2(s,i) 

The diagonal elements of hki{s, t) are of the form 

-4s -2t + l2Ml 



(76) 



hoo{s,t) 
hii{s,t) 



37rs(s - 4M2)(s -t){s + t~ AM^) 

-3s - ^ + 8Afg 

27rs(s - 4M2)(s -t){s + t~ AM^) 
-7s~5t + 2AMl 



, (77) 



67rs(s - 4M2)(s ~t){s + t- AM^) 
The off-diagonal ones are given by 

Cke 



'^"''-''^^ 6^s(s-4M2)(s + i-4M2) 

coi = 6, co2 = -10, 

cio = 2, ci2 = -5, 

C20 = -2, C21 = -3. 



, (78) 



The kernels in the second line of (|75p can be expressed in 
terms of those listed above: 

^'^'''^ = i^AMl ' 

Co = 2 , Ci = 1 , C2 = -1 . 

As in the case of S{t), the integrands in ([75)) display 
a fictitious double pole for s = AM^ — t, but the residue 
of the singularity vanishes for a crossing symmetric in- 
tegrand. At i = 0, the three sum rules (TM)) reduce to a 
single one, in fact one that had been made use of already 
earlier F°l 

Since the contribution from the S- and P- waves is man- 
ifestly crossing symmetric, these waves do not show up in 
the above sum rules. For the S-waves, which generate a 
independent contribution, that is manifest: these waves 
drop out, in ImT(s,t), in ImT^(s,0), as well as in the 
terms involving the partial derivative with respect to t. 
The P^-wave, on the other hand does generate a contri- 
bution to these quantities - it drops out only in the sum, 
on account of the relation (|75|) between the kernels. 



D Details of minimization procedure 

In the present appendix, we specify the discrepancy func- 
tion that enters the procedure used to solve the sum rules 
and duality conditions, as described in section O For a 



given input, these amount to constraints imposed on the 
Regge output variables, 

^out = {pi,bp,cp,(3p{0),ri,bp,Cp,l3e{0),be,Ce}. (80) 

We evaluate the ^-dependent conditions on a finite number 
of points, using the lattice 



t„ = nt-i 



-0.1GeV^ 71 = 0,. 



,6. 



(81) 



The sum rules and duality conditions then amount to a 
finite number of relations relations of the form 

^PW + 5R(„Rj^0, r = l,2,... 

The part represents the contribution from the region 
where we use the partial wave representation, while the 
remainder is calculated with the Regge parametrization. 

The duality conditions for the total cross sections are 
solved as algebraic constraints for the variables pi , r i , /3e (0) . 
The remaining variables in v^^^ are determined by mini- 
mizing the discrepancy function 



(82) 



See appendix B.2 in equation (B.7). 



where the sum extends over all constraints except the du- 
ality conditions at t = 0. The factors specify the weights 
of the remaining conditions in the minimization process. 
We determine these with the response of the term to 
variations of the partial wave input: vary the input param- 
eters within the estimated uncertainty range, calculate the 
responses generated by these variations in 5^^, evaluate 
the sum Ar over the squares of these responses and set 
\r = l/Ar- The choice of the weights is not of crucial 
importance, because the various sum rules and duality 
conditions are consistent with one another. The particu- 
lar choice made homogenizes the set of relations, giving 
little weight to those constraints that are sensitive to the 
uncertainties in our partial wave representation and thus 
contain little information about the Regge output param- 
eters, which we determine by minimizing X^(^(^t)- 
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